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NOTATION 


Area  of  a  lifting  surface 

Area  of  the  helical  surface  behind  the  trailing  edge  of  a 
lifting  surface  extended  to  infinit\ 

Area  between  a  superimposed  lifting  line  and  the  trailing 
edge 

Coefficients  as  defined  in  Equation  ll7] 

Coefficients  as  defined  in  Equation  l2ll 
Ordinate  of  mean  line 
Diameter  of  propeller 

Advance  coefficient  based  on  ship  velocity 

Chordwise  coordinate  of  any  point  with  respect  to  reference 
line 

Chordwise  coordinate  of  leading  edge 
Cnordwise  coordinate  of  trailing  edge 
Chord  length 

Revolutions  per  unit  time 

Point  on  lifting  surface  at  which  mean  line  is  sought 
Pitch-diameter  ratio 

Distance  vector  between  two  points 
Radial  coordinate 
Radial  coordinate  of  hub 

Radial  coordinate  of  point  F 

Length  vector  of^an  elemenlarv  bound  vortex  line  on  a 
lifting  surface;  clr/\(ir\  is  a  unit  vector  tangent  to  the 
bound  vortex  line 

Length  vector  of  an  elementary  free  cortex  line,  (is/\rd(^\ 
is  a  unit  vector  tangent  to  the  free  vortex  line 


Nondimensional  axial  and  tangential  componenU;,  respec- 
ti\ely,  of  induced  velocity  from  li^iing-line  calculation 
based  on  ship  velocity 

Nondimensional  inflow  \elocity  based  on  ship  speed 

Induced  velocity  vector  at  point  P  due  to  radial  bound 
circulation 

Axial  and  tangential  components  of 

Induced  velocity  vector  at  point  P  due  to  free  circulation 

Axial  and  tangential  components  of  Vp 

Velocity  vector  at  a  lifting  line  induced  by  trailing  free 
vortex  sheets 

Total  induced  velocity  vector  at  point  P 

Total  induced  velocity  vector  at  point  P  relative  to 

Ship  velocity 

Nondimensional  local  longitudinal  velocity  based  on  ship 
velocity 

Nondimensional  induced  velocity 

Nondimensional  axial  and  tangential  components,  respec¬ 
tively,  of  induced  velocity  v  based  on  ship  velocity 

Nondimensional  induced  downwash  at  point  P  normal  to 
helical  chord  line  based  on  ship  velocity 

Axial  and  tangential  components  of  Vg  within  a  strip 
*  between  and 

Axial  and  tangential  components  of  Vp  within  a  strip 
between  Tj  and 

Cartesian  coordinates 
Chordwise  coordinate 
Local  wake  fraction 

Radial  coordinate  as  defined  in  Equation  [17] 

Number  of  propeller  blades 


/3 

Advance  angle 

Hydrodynamic  pitch  angle 

r 

Nondimensional  circulation  =  circulation  ~DV^ 

r(r) 

Nondimensional  bound  circulation  distribution  along  a 
lifting  line 

Nondimensional  free  circulation  distribution  on  free  vortex 
sheet  trailing  a  lifting  line  or  a  lifting  surface 

Vfir,  6) 

Nondimensional  free  circulation  distribution 

r,(r,  0) 

Nondimensional  bound  circulation  distribution  on  a 
lifting  surface 

r^Cr,  6) 

Nondimensional  free  circulation  distribution  on  a 
lifting  surface 

y 

Circulation  per  unit  angle  =  — ~ 

Oj 

Half  width  of  a  strip  between  and 

X-Xq 

e 

Angular  coordinate  =  L  cos 

Angular  coordinate  of  leading  edge 

Angular  coordinate  of  point  P 

Angular  coordinate  of  trailing  edge 

Oj 

Total  chord  length  in  terms  of  angular  coordinate 
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ABSTRACT 


This  report  presents  a  method  of  the  propeller  camber  calculations  for 
the  case  of  uniform  chordwise  load  distribution  employing  the  lifting-surface 
theory.  This  is  essentially  a  refinement  of  Pien’s  published  work  on  propeller 
lifting-surface  theory.  Mathematical  development  |>ertinent  to  the  numerical 
computation  is  reviewed.  Details  of  the  computational  procedure  and  method 
are  outlined,  and  results  of  sample  calculations  are  included.  A  detailed 
instruction  for  preparation  of  input  data  for  the  computer,  samples  of  the 
computer  input  data  and  output,  and  the  FORTR.AN  listing  of  the  computer 
program  are  also  given. 


INTRODUCTION 

In  Reference  1,  Pien  proposed  an  improved  method  for  the  calculations  of  marine  pro¬ 
pellers  by  applying  the  lifting-surface  theory.^  Included  in  Reference  1  was  an  outline  for 
the  numerical  computation  by  means  of  a  high-speed  computer,  IBM  709.  The  computation 
was  programmed  for  the  case  of  uniform  chordwise  load  distribution.  After  his  paper  was 
published,  a  model  propeller  (Model  3916)  was  designed  by  using  this  improved  method  and 
was  tested;  the  results  showed  a  thrust  deficiency  of  approximately  7  percent. 

In  the  fall  of  1962  an  iCtive  propeller  research  project  was  planned  and  initiated  by 
Pien.  The  overall  progran  included  the  formulations  and  the  numerical  evaluations  for  the 
following  cases: 

1.  Propeller  camber  calculation  method  for  uniform  chordwise  load  distribution, 

2.  Propeller  cainber  calculation  method  for  arbitrary  (nonuniform)  chordwise  load 
distribution, 

3.  Prediction  of  performance  for  a  given  propeller  in  steady  inflow,  and 

4.  Prediction  of  performance  for  a  given  propeller  in  unsteady  inflow. 

This  report  deals  with  the  first  case,  the  camber  calculation  method  for  uniform  chord- 
wise  load  distribution.  With  the  overall  project  in  view,  the  numerical  computation  procedure 
and  the  programing,  reported  in  Reference  1,  were  revised.  This  new  program  incorporates 
features  that  will  provide  a  greater  degree  of  flexibility  to  facilitate  handling  of  a  larger 
number  of  parameters  and  to  serve  as  a  basic  program  for  the  subsequent  works. 

A  brief  discussion  of  the  formulation  of  the  problem  and  a  detailed  description  of  the 
numerical  computation  procedure  arc  given  in  this  report.  -  Results  of  a  number  of  examples 
with  different  design  parameters  such  as  number  of  blades  and  area  ratios  are  also  included. 


References  are  listed  on  page  61. 


A  detailed  instruction  for  preparation  of  input  data,  and  the  FORTRAN  listing  of  the  computer 
program,  as  well  as  samples  of  computer  input  and  output  data  are  given  in  the  appendixes. 

FORMULATION 

The  detailed  formulation  of  the  problem  is  given  in  Reference  1.  For  the  sake  of  com¬ 
pleteness  of  presentation,  a  brief  discussion  of  the  mathematical  development  of  the  problem 
and  the  derivation  of  the  pertinent  equations,  on  which  the  numerical  computations  are  based, 
are  given. 

COORDINATE  SYSTEM 

A  Cartesian  (a;,  y,  2)  coordinate  system,  shown  in  Figure  1,  is  fixed  on  the  propeller. 
The  z-axis  is  the  axis  of  the  propeller  shaft  with  its  positive  pointing  upstream,  the  y-axis 
is  an  axis  on  the  first  blade  arbitrarily  selected  so  as  to  pass  through  the  midchord  at 
r  =  0.7R  with  its  positive  outward,  and  the  2-axis  is  the  third  axis  and  has  its  positive  in 
accordance  with  the  right-hand  rule.  Another  coordinate  system  is  also  shown  in  Figure  1. 

It  is  a  cylindrical  coordinate  system  in  which  a  point  in  space  can  be  identified  by  (r,  0,  x) 
where  r  is  the  radial  coordinate,  6  is  the  angular  coordinate  (positive  clockwise  looking  for¬ 
ward),  and  X  is  the  same  as  defined  above.  Along  a  helical  line  the  axial  coordinate  can  be 
expressed  in  terms  of  pitch-diameter  ratio  p(r),  i.e.,  x  =  p(r)d/n  as  shown  in  Figure  1.  This 
cylindrical  coordinate  system  is  used  in  the  fwmulation  of  the  basic  equations. 

For  convenience  in  the  numerical  computations  a  “moving”  cylindrical  coordinate  is 
devised.  This  system  has  r,  and  ^  as  its  coordinates,  where  r  is  the  radial  coordinate, 
ff>  is  the  new  angular  coordinate  and  is  defined  as  <f>  =  6  -  9^,  and  f  is  a  new  axial  coordinate 
and  is  defined  as  ®  -  z^,  where  (r,  0,  z)  and  (r^,  9^  ,  z^)  identify  a  field  point  Q{r,  0,  z) 
and  a  control  point  0^,  z^),  respectively,  in  the  “stationary”  cylindrical  coordinate 

system.  Also  defined  are  the  unit  vectors,  i,  /,  k,  in  the  directions  of  r^,  0,  and 
respectively. 

In  this  new  “movi^Jg”  cylindrical  coordinate  system  the  control  point  F(rp,  0^,  z^)  is 
identified  by  the  vector  i  r^,  and  the  distance  vector  R  between  a  field  point  Q(r,  0,  z)  and  P 
is  expressed  as  follows: 

R  =  i  -  r  cos  4>)  -  j  r  sin  0  -  kff>p  {r)/Tt 

Similarly,  the  vectorial  expressions  for  the  element  of  length  of  the  bound  vortex  line  on  the 
lifting  surface,  dr,  and  fc»-  the  element  of  length  on  the  free  vortex  line,  ds,  are  as  follows: 

dr  =  dr  (i  cos  ^  sin  <f>  +  k  if)  p  '{ry/n) 
ds  =  d(f){-ir  sin  <f)  +  jr  cos  if)  +  kp{r)/n) 

These  two  expressions  will  be  used  in  the  evaluation  of  the  induced  velocities. 
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BASIC  EQUATIONS 


It  is  assumed  that  the  propeller  is  operating  in  incomfX’essible,  inviscid  fluid.  The 
thickness  of  the  propeller  blade  is  assumed  to  be  zero.  The  theory  postulates  that  the  pro¬ 
peller  blade  may  be  represented  by  a  lifting  surface  of  the  same  outline.  On  this  surfac*', 
vorticity  distribution  is  so  placed  and  with  such  intensity  that  for  a  given  design  condition, 
a  certain  force  distribution  is  developed. 

Let  r,  d,  and  x  be  the  coordinates  of  a  field  point  Q  (see  Figures  1  and  2)  where  the 
intensity  of  bound  vorticity  is  denoted  as  (r,  6).  Then  the  following  relationship  holds: 

Ji(n 

r^(r,d)dd  =  r(r)  [1] 

(r) 


where  chordwise  coordinates  d^(r)  and  0^(r)  define,  respectively,  the  leading  and  trailing  edges 
of  the  lifting  surface,  and  F  (f)is  the  radial  bound  circulation  distribution  which  can  be  com¬ 
puted  from  lifting  line  theory  to  satisfy  certain  design  thrust  or  torque  conditions. 

Because  of  the  variation  of  F (r)  along  the  radial  direction,  a  free  vortex  sheet  is  shed 
and  trails  behind.  It  is  assumed  that  this  trailing  vortex  sheet  is  a  general  helical  surface 
formed  by  helical  lines,  neglecting  the  radial  contraction  of  the  slipstream.  As  an  approxi¬ 
mation  the  helical  surface  is  extended  to  the  leading  edge  of  the  lifting  surface  and  the  lifting 
surface  is  assumed  to  be  a  part  of  this  helical  surface.  In  the  lifting-surface  representation 
there  exist  two  vortex  systems:  the  radial  bound  system  on  the  surface  with  an  intensity  of 
r'r  system  with  an  intensity  of  F^  (r,  (?).  The  latter  may  be  considered  as 

consisting  of  two  parts:  the  part  that  is  on  the  lifting  surface  with  an  intensity  designated 
as  F^  (f,  0),  and  the  part  that  trails  behind  the  trailing  edge  in  the  slipstream  which  has  an 
intensity  of  F^  (r).  In  the  slipstream,  the  circulation  F^  (r)  is 


r,  =  r(r)dr 

dr 

Combining  Equations  [1]  and  [2],  we  have 


rc(r)  =  - 


r^{r,e)de 


dr 


[2] 


[31 
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The  free  circulation  on  the  lifting  surface  may  be  expressed  as 


$1 

(r,  0)  =  -  ^  j  r,  (r,  6)  ddj 


The  induced  velocity  at  any  point  P  on  the  lifting  surface  due  to  these  two  vortex  sys¬ 
tems  is  derived  by  the  use  of  the  law  of  Biot-Savart.  The  resulting  general  expression  is  as 
follows: 


6)  dO 


drx.  R 


ds  y-R 

6)  -  [4] 


where  Vp  =  total  induced  velocity  vector  at  point  P  on  the  lifting  surface, 

R  =  distance  vector  between  a  point  P  and  a  point  where  elementary  vortex  is  located, 

ds  =  length  vector  of  an  elementary  vortex  line  in  the  free  system, 

dr  =  length  vector  of  an  elementary  bound  vortex  line  on  the  lifting  surface, 

=  area  of  the  lifting  surface,  and 

^2  =  area  of  the  helical  surface  behind  trailing  edge  of  the  lifting  surface  extended  to 
infinity  downstream. 

The  second  term  of  Equation  [4]  may  be  expressed  as  follows: 


ill 


I7(r,  0) 


ds  X  R 


■ill 


(r,  0) 


dsx  R 


•ill 


ds  X  R  ,  , 
r.w — —  15] 


+  A2 


Substituting  Equation  [5]  into  Equation  [4]  and  adding  and  subtracting  a  quantity  of 


—  ff  -  - 

4.7  JJ  dr 


ds  X  R  ,  j 

r  (r)  dr  • - ,  where  is  the  enclosed  area  between  the  line  drawn  on 


the  lifting  surface  through  the  centerline  of  shaft  and  the  point  P  and  the  trailing  edge,  we 
obtain  the  following  expression  for  the  induced  velocity: 


1  r  ff 

L  J  J  A 


0)  d0 


dr  X  R 


r^Cr,  0) 


11,  ff  ± 

JJ,  dr 


r{r)dr 


ds  X  R 


r  r 
in  J « 


d  ds  X  R 

—  r(r)(fr  - —  + 

dr  1^3 


d  ds  X  R 

-  -  r{r)dr  - — 
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t? 


since  P  (r)  =  -  —  F  (r)  dr  by  Equation  [2]. 

dr  _ 

The  second  square  bracketed  term  is  a  function  of  r  only  and  is  denoted  as  .  It 
expresses  the  induced  velocity  at  a  lifting  line  with  a  bound  circulation  F  (r).  This  term 
plays  a  very  important  role  in  our  problem  of  hydrodjnamic  aspect  of  propeller  design  using 
lifting  surface  theory.  Consider  the  case  in  which  the  pitch  of  the  lifting  surface  is  constant. 
If  a  lifting  line  with  bound  circulation  F  (r)  is  superimposed  on  the  surface  through  the  control 

point  P,  then  it  w'ill  shed  a  vortex  system  with  a  circulation  of  -  —  F(r)rfr  and  result  in  an 

_  dr  _ 

induced  velocity  at  P  which  is  identical  to  The  values  of  this  induced  velocity  at 
various  radii  may  be  computed  using  lifting  line  theory  based  on  a  chosen  radial  load  distri¬ 
bution.  From  these  values  the  helical  surface  near  the  {wopeller  disk  is  then  determined, 
upon  which  we  place  the  lifting  surface  with  a  chosen  blade  outline.  Since  is  independent 
of  0,  it  is  convenient  to  use  it  as  the  reference  velocity  for  the  induced  velocity  due  to  a 
lifting  surface.  Let  us  define  the  total  induced  velocity  at  point  P  relative  to  Fj^,  Fp^  = 

yp-\^ 

If  Fp^  is  zero  everywhere,  then  Fp  =  F^  and  the  streamline  on  the  lifting  surface  will 
coincide  with  the  helical  chord  at  each  radius.  If  Fp^  is  not  zero  everywhere,  then  it  will 
induce  a  mean  line  with  respect  to  the  helical  chord  at  each  radius.  Since  we  are  only  inter¬ 
ested  in  the  mean  line  with  respect  to  the  helical  chord,  is  the  only  item  that  needs  to  be 
computed  and  we  can  disregard  Fp.  Carrying  this  idea  through,  we  obtain  from  Equation  [6] 
the  following  expression  for  Fp^: 


—  1 
F„ ,  =  — 


PL 


in 


C  C  dr  X  R  f  f  ds  X  R 

J J  F  (r,  d)  dd  +  J  J  F0(r,  6) 


II 


R' 


ds  X  R 


—  F(r)rfr 
dr  1^3 


=  Fp  +  Fp 


[7] 


where 


Fo  = 


in 


F,(r,  d)de 


dr  X  R 


[8] 
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If  we  compare  Equation  [7]  with  Equation  [6],  it  is  obvious  that  the  amount  of  numerical  com¬ 
putation  is  reduced  quite  appreciably.  In  other  words,  our  problem  has  been  greatly  simplified 
by  taking  advantage  of  these  shortcuts.  This  is  the  unique  feature  of  Pien’s  approach  to  the 
problem. 

It  should  be  mentioned  here  that  the  foregoing  procedure  is  true  only  for  the  case  where 
the  pitch  of  the  lifting  surface  is  constant.  If  the  pitch  changes  with  radius,  the  free  vortex 
sheet  of  a  lifting  line  may  not  exactly  coincide  with  that  of  the  lifting  surface.  However,  the 
influence  of  this  on  the  induced  downwash  is  rather  small,  and  this  procedure,  therefore,  may 
also  be  applied  to  the  cases  where  the  pitch  is  not  constant. 

The  integrals  on  the  right-hand  side  of  Equation  [7]  now  have  bounded  limits  of  inte¬ 
gration  and  can  be  handled  conveniently. 

From  the  axial  and  tangential  components  of  the  relative  induced  velocity  the 
downwash  distribution  ■w^,  normal  to  the  chord,  is  obtained.  By  integrating  this  downwash 
distribution  curve,  the  required  mean  line  relative  to  the  chord  is  obtained. 


^0^ 


2V  cos  ^ 


■  j  '-.(v 

»  ''a 


[10] 


EQUATIONS  FOR  NUMERICAL  COMPUTATION 

The  major  portion  of  the  work  is  in  the  evaluation  of  the  induced  velocities.  By  using 
Equation  [7]  as  a  basis,  the  following  working  equations  for  numerical  computations  are 
develop>ed. 


Induced  Velocity  at  Point  P  Due  to 
Bound  Vortex  System,  Vg 

The  induced  velocity  at  any  control  point  P  due  to  the  bound  vortex  system  is  expressed 
by  Equation  [8].  To  reduce  it  into  a  working  form,  we  use  the  “moving”  cylindrical  coordinate 
system  for  convenience.  Making  use  of  the  expressions  for  dr  and  R  developed  previously,  we 
obtain  the  following  equation  for  induced  velocity  due  to  the  bound  vortex  system: 


r,(r,^)  _ 

-  {  i  {-4>  sin  <f>p{T)/:T  +  r<^  sin  <f>p'{r)/n] 

+  /  cos  <f>p{T)/77  +  4>(rQ-f  cos  4>)p'(r)/7T] 

+  kl-r^  sin  <^]\  d<f>dr  [11] 

Induced  Velocity  at  Point  P  Due  to 
Free  Vortex  System,  Vp. 

In  a  similar  manner  we  put  Equation  [9]  into  a  working  form  by  substituting  the  expres¬ 
sions  for  R  and  ds  developed  previously.  It  follows  that 


where  R  =  [r^  -  2 r^r  cos 


1  - 

—  cos  <5  +  sin  <^)p(r)/7i 

R^ 


*1  — 


+  i  (fg  -  f  COS  <fi  -  T(f>  sin  ^)p(r)/n- 
+  k  r  {t  -  Tq  cos  ^)]  [r ^(r,  0)  +  P  '(r)  rfr]  d<f> 


[121 


where  r'(f)  =  —  r(r)  equals  zero  when  ^  >  0. 

dr  _  _ 

The  total  induced  velocity  at  point  P  relative  to  is  therefore  the  sum  of 

Equations  [11]  and  [12].  For  the  present  problem  the  components  that  are  to  be  evaluated 
are  the  axial  (A)  and  the  tangential  (;)  components. 

To  reduce  Equation  [12]  it  is  necessary  to  obtain  an  expression  for  the  intensity  of 
the  free  vortex  system  that  is  on  the  lifting  surface,  r^(r,  ^).  This  may  be  obtained  by 
performing  the  differentiation  in  Equation  [3]. 


r  fr)  =  -rfr 


>r) 


—  r^(r,  e)'dd  -  r^(r,  e^)  ^  0p(r)</r  +  r,(r,  0;)  0,(r)  dr 


Then,  within  the  lifting  surface,  the  free  circulation  can  be  expressed  as  follows: 


0)  =  -  dr 


0,(0 


—  r,{r,  0)  d0  -  r^(r,  0)  —  0^(r)  dr  [13] 
dr  or  ‘ 
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For  the  case  of  uniform  chordwise  load  distribution,  the  bound  circulation  is  uniformly  distri¬ 
buted  along  the  chord.  Its  intensity  per  unit  angle  is  obtained  from  Equation  [1]: 


<f>)  = 


r(r) 

Oi(r)-d^{r) 


ISl 

Ojir) 


[14] 


Then,  Equation  [13]  can  be  further  simplifiod  and  becomes 

4>)  =  -  *  W  -  ^)  —  y(r)  +  y(r)  j  [15l 

With  Equation  [13]  substituted  into  Equation  [11],  and  Equation  [15]  into  Equation  [12], 
we  obtain  the  following  expressions  for  the  axial  and  tangential  components  of  the  induced 
velocities; 


sin  ^ 

“  '■fl  yW  -  d<l>dr 


[16a] 


1  rr  1  <f> 

^  7^  ~7  ^  d<f>dr 

d  J  A  R^ 


[16b] 


''■■iil- 


1  /Y  '■^'■"^0  *^) 

- ; -  [[(^iW-0)  y'(r)  +  y(r)0/{r)]  -T'{^\di>dT 


[16c] 


Pif) 


(fg-r  cos  sin  0)  }[(<^,(r)-0)y'(r)  +  y(r)  0f(r)] 

-r'(r)U<^rfr 

[16d] 


for  0  >  0  the  r'(f)  term  in  Equations  [16c]  and  [16d]  is  equal  to  zero. 

The  total  axial  component  of  the  induced  velocity  is,  therefore,  the  sum  of  Equations 
[16a]  and  [16c],  and  the  total  tangential  component  is  the  sum  of  Equations  [16b]  and  [16d]. 
These  equations  are  the  basic  working  equations  that  are  to  be  evaluated  numerically. 

The  Equations  [16a]  through  [16d]  are  developed  for  the  induced  velocity  at  a  point  P 
on  the  first  blade  due  to  a  single  blade,  the  first  blade  itself.  Since  all  the  blades  may  be 
represented  by  lifting  surfaces,  there  will  be  contributions  due  to  each  of  these  surfaces  to 
this  induced  velocity.  The  coordinates  of  a  field  point  on  m*’’  blade  corresponding  to  a  field 
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point  (r,  on  the  first  blade  at  the  same  radius  is  (r,  where  =  +  2n-(jn-l)/Z, 

m  =  1,  2,  3,  . . .  ,  Z,  assuming  the  blades  are  arranged  symmetrically  with  equal  angular 
spacings.  To  derive  an  expression  for  the  induced  velocity  due  to  the  blade,  it  is  only 
necessary  then  to  replace  sin  ^  by  sin  and  cos  <f>  by  cos  in  Equation  [16]  since  the 
blades  are  assumed  to  be  identical  to  each  other.  Thus,  for  example,  the  axial  component  of 
the  induced  velocity  due  to  the  bound  vortex  system  of  the  blade  is  derived  from  Equa¬ 
tion  [16a]: 


1 

An 

It  follows  that  if 
all  the  blades: 


If 


Tq  <r) 


sin  <f> 


d(f)dr 


we  sum  up,  we  would  obtain  an  expression  for  this  induced  velocity  due  to 


1  ^  r  r  sin  d> 

=  i  I  //-'oHO- - ^ 

m  =  1  [fp  -  2rQ  r  cos  <f>^+  ^ 

Similarly,  the  various  other  induced  velocity  components  due  to  all  the  blades  may  be  derived 
for  the  three  remaining  equations. 


NUMERICAL  COMPUTATION 

As  stated  previously,  the  radial  bound  circulation  distribution  and  the  hydrodynamic 
pitch  ratio  distribution  are  obtained  from  the  lifting-line  theory  based  on  certain  design 
conditions.  (At  the  David  Taylor  Model  Basin  the  program  entitled  “AML  Problem  No.  106  - 
Moderately  Loaded  Propeller  Using  Induction  Factors”  may  be  used  for  this  purpose.)  These 
quantities  along  with  the  blade  outline  and  the  chord-length  ratios  L^/D  and  Lj/D,  which  are 
chosen  based  on  the  considerations  of  cavitation,  blade  element  efficiency,  and  strength,  are 
used  as  input  at  given  radii  for  the  numerical  computations. 

For  the  purpose  of  interpolating  these  quantities  at  radii  other  than  those  specified  in 
the  input,  trigonometric  series  representations  in  finite  number  of  terms  of  these  quantities 
are  made  by  using  the  expression 

,  AT 

A  =  ^  cos  n  y  [17] 

n  =  0 
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where  X  denotes  any  of  the  preceding  quantities 
y  =  radial  coordinate  defined: 

cos  y  =  (1  +  -  2r)/(l  -  r^) 

=  coefficients  determined  from  input  data 
N  =  number  of  terms. 

From  this  expression,  the  quantity  X  can  be  calculated  at  any  desired  radial  position.  The 
derivatives  of  these  quantities  with  respect  to  r  can  be  easily  expressed  as 


—  X  2 

dr  (1  -  f^)  sin  y 


N 

^  71  sin  71  y 

n  =  0 


[183 


and  they  can  also  be  computed  at  any  desired  radius. 

For  each  radius  at  which  the  blade  element  mean  line  relative  to  the  chord  line  is 
sought,  various  values  of  6^  are  chosen  so  that  the  distributions  of  components  of  the  induced 
velocity  can  be  obtained.  At  any  point  P  (r^,  6^)  the  induced  velocities  are  computed  according 
to  Equations  [I6a3,  [16b],  [16c],  and  [i6d]. 

It  should  be  noted  that  in  these  equations  the  integrand  has  a  singular  point  at  r  = 
and  6  =  6q{4>  -  0).  The  integrations  are  carried  out  by  isolating  the  singular  point.  The 
integrals  are  evaluated  numerically  in  the  regions  of  |r  -  >  5^,  where  is  a  quantity  that 

can  be  made  as  small  as  practicable.  Within  the  strip  \t  -  f^l  <  5^,  the  region  of  integration 
is  divided  into  two  parts.  The  first  part  is  over  the  portion  |0|  >  <o’  second  part  is 

over  the  portion  1^1  <  <g.  The  first  portion  again  can  be  carried  out  numerically.  The  second 
portion,  |r  -  r^l  <  8^  and  |«^|  <  «q,  is  not  evaluated,  for  when  5^  and  4^  are  very  small,  the 
contribution  to  the  camber  due  to  this  portion  is  negligible  even  though  the  induced  velocity 
may  be  infinite  at  some  of  the  singular  points.  (See  Appendix  A.) 

The  numerical  computation  for  the  narrow  strip  of  the  lifting  surface  within  which  the 
control  point  P  is  located  is  handled  differently  from  the  remaining  portions  of  the  lifting 
surface.  Thus,  to  facilitate  the  numerical  calculation  the  lifting  surface  A  is  divided 
geometrically  into  three  zones  (Figure  2):  the  “mid-zone,”  which  is  the  narrow  strip  with  a 
width  of  28 Q  from  r  -  5^  to  the  “tip-zone,”  which  is  the  area  between  and  the 

tip;  and  the  “hub-zone,”  which  is  the  area  between  the  hub  and  r  -  8^.  The  division  of  zones 
varies  with  the  radial  position  of  the  point  P.  The  computation  procedures  in  the  tip- zone 
and  the  hub-zone  are  identical.  For  numerical  integration  the  Simpson  method  is  used.  In  the 
tip-  and  hub-zones,  the  chordwise  integration  is  performed  first,  and  in  the  narrow  mid-zone, 

which  has  a  width  of  25^,  the  radial  direction  integration  is  done  first.  Within  this  mid-zone 

d 


the  pitch  of  the  chord  p(r),  —  p(r)  and  the  circulation  functions  F  (r,  ^),  rfl(f,  0),  and 
d  , 

—  r(r)  are  considered  to  be  independent  of  the  radial  coordinate,  and  the  mean  values 
dr 
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evaluated  at  the  midwidth,  r^,  of  the  strip  are  used.  Then  the  integrations  in  Equations  [16aj, 
[16b],  [16c],  and  [16d]  can  be  carried  out  for  the  interval  between  and  +  5^: 


} 


ilrj 


(r  -  Tq  cos  <f>)  sin  <f> 

[r^  sin^  0  +  (<f>p{rQ)/7Tf]R 


d<f> 


[19a] 


^0-^0 


•^t(0 


0 


y(^o) 


r^s\n^<l>+(4,pirQ)/ny 


r-Tp  CCS  4> 

(p(ro)  COS  0  +  ro  p  '(r^) - - - 


+  p' cos  <!> 


^0  +  (<f>p(roV^f  -  '■'■q  <f> 


R 


dtf) 


’’o  ^0 


[19b] 


W  =  i  j 


In  (r  -  r^cos  <f)  +  R)  -  - 


d<b 


[19c] 


f 

^TT  J 


-iSF) 


p(ro)  j  rr^  (sin^  ^  ^  sin  <p  cos  +  sin  ^  P(^q)/”’)^  (cos  0  +  ^  sin 


[r^  sin^  <?!.  +  {<l>p{r^)/n)-]R 


'■o  +  ^o 
)>-  d<i> 
^0  ~  ^0 


[19d] 


where  SF  =  -  <f>)  y'i^o)  +  y('’o)  ^/('’o)  ~  ^ '(’’o) 

r'term  equals  zero  for  dt  >  0.  These  chordwise  integrations  are  then  computed  numerically 

by  the  Simpson  method. 

The  contribution  of  other  blades  is  computed  in  the  sanre  manner  with  the  blade  repre¬ 
sented  by  a  lifting  surface.  In  the  computation  for  the  m**’  blade  the  sin0  and  cos  0  terms 
are  replaced  by  sin  (0  +  2n{m  -  1)/Z)  and  cos  (0  f  2rr(m  -  1)/Z)  respectively,  where  Z  is  the 
number  of  blades,  m  =  2,  3,  . . .  ,  Z. 
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The  total  axial  and  tangential  induced  velocity  components  and  are  then 
expressed  as 


=  '^ha  +  ^fa  ^  ^fa 

+ 

where  the  primes  denote  the  velocity  components  due  to  other  blades. 

From  these  induced  velocity  components  the  downwash  distributions  relative  to  the 
chord  is  calculated  by  using 


('■q’  ('“o*  ^o)  /3,.  -  v^(r^,  6^)  sin  /S. 


[20] 


For  calculating  the  blade  mean  line,  the  downwash  is  expressed  in  a  power  series  in 

terms  of  a  new  nondimensional  angular  coordinate  —  as  follows: 

S 


M 


[21] 


n  =  0 


where  the  integer  M  equals  the  number  of  chordwise  coordinates  minus  one,  and  the  coeffi¬ 
cients  (r^j)  are  evaluated  from  (Af  +  1)  simultaneous  equations.  With  this  scheme  the  down- 
wash  at  any  point  can  be  obtained  by  either  interpolation  or  extrapolation,  and  the  integration 
of  Equation  [10]  can  easily  be  performed  to  obtain  the  blade  mean  line  ordinates: 


C 


m 


i  '  0'  ^  n 

V{r^)  \  ~1 

n  =  0 


[22] 


where  V  {r^)  is  the  mean/esultant  velocity.  The  mean  resultant  Velocity  is  computed  based 
on  the  following  geometric  relationship;  see  Figure  3. 

where  =  nondimensional  local  longitudinal  velocity  based  on  ship  velocity; 

“a’  ~  nondimensional  axial  and  tangential  components,  respectively,  of  induced 

velocity  from  lifting-line  calculation  based  on  ship  velocity;  and 

=  advance  coefficient  based  on  ship  velocity. 
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THE  COMPUTER  PROGRAM 


The  numerical  computation  has  been  programed  into  the  IBM  7090  machine.  The  com¬ 
puter  program  is  designated  as  Applied  Mathematics  Laboratory  Problem  XPLU.  First,  the 
program  calculates  the  angular  positions,  0,  of  field  points  relative  to  any  point  P  at  which 
the  blade  element  mean  line  is  sought  according  to  the  following  relationships: 

n  ^  ^l,n-k+l  ~  ^l.k 

with 

and 

^k.  1  =  ® 

where  k  =  1,  2,  3,  . . .  and  n  =  2,  3,  4,  ...  are  the  indices  for  the  radial  and  chordwise  stations, 
respectively. 

(A0)^.^  is  the  minimum  spacing  specified,  and  (n  -  1)  (A^)  is  spacing  between  the 
and  the  (n  -  I)***  stations.  When  it  exceeds  a  prescribed  maximum  value  in  the 

input,  it  will  take  on  the  value  of  i.e-,  the  spacing  A<;6  is  limited  to  <  A^ 

<  i^4>)aiax'  shows  that  the  field  points  are  so  located  that  they  have  smaller  spacings 
in  the  close  neighborhood  of  point  P,  and  the  spacing  increases  as  the  distance  between  P 
and  the  field  point  increases.  The  fineness  of  the  network  varies  depending  upon  the 
and  specified.  For  each  of  these  angular  positions  the  sine  and  cosine  of  the  angles 

are  computed,  compiled  into  tables,  and  stored.  These  tables  are  used  for  providing  ready 
information  of  0,  s\n<p,  and  cos  <f),  required  in  the  later  stage  of  the  computation.  Since  only 
one  set  of  angles  is  computed,  a  common  pattern  of  network  of  field  points  for  any  point  P  is 
thus  established.  The  specified  maximum  range  of  angles  for  which  the  tables  are  made  is 
sufficiently  large  so  that  they  may  be  used  for  any  propeller  of  the  conventional  configuration. 
Next,  the  program  calculates  the  coefficients  in  Equation  [17]  with  the  input  data  supplied. 
The  number  of  terms,  N,  depends  on  the  number  of  intervals  of  radii  r  in  the  input  data  which 
can  be  chosen  arbitrarily.  A  9-interval  of  r  is  recommended,  and  the  present  program  is 
programed  for  9  intervals. 

The  computer  then  selects,  according  to  the  input,  the  radial  coordinate  (starting  from 
the  innermost  radius)  of  point  P  at  which  the  blade  element  mean  line  is  sought.  The 
radial  stations  where  the  elemental  vortex  is  located  in  each  of  the  tip-  and  hub-zones  are 
calculated  according  to  the  spacings  specified  in  the  ipput  data. 

At  each  of  those  radiai  stations  r,  quantities  that  are  functions  of  rand  the  derivatives 
required  in  the  computation  of  the  induced  velocities  are  then  calculated  and  stored. 

For  a  given  radial  station  r  the  factors,  which  are  functions  of  rand  (f>,  appearing  in 
the  expressions  of  the  induced  velocities  are  first  calculated  for  various  angles  at  r  and  then 
are  stored. 
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The  computer  proceeds  to  select  a  chordwise  coordinate  (starting  from  the  leading 
edge)  0Q  of  point  P  according  to  the  specified  intervals.  The  chordwise  stations  ^  of  the 
field  points  are  then  computed  with  the  aid  of  the  table  of  ^  compiled  at  the  beginning  of  the 
program  at  the  corresponding  radial  station.  The  number  of  chordwise  stations  depends  on  the 
difference  between  r  and  and  the  blade  outline  of  the  particular  propeller  involved.  This 
may  be  better  explained  by  a  gfraphic  illustration.  Figure  4  shows  a  pattern  of  a  network  of 
field  points  for  point  P,  whe-s  both  the  radial  and  the  chordwise  stations  are  symmetrical 
about  the  point  P.  In  Figure  5,  a  propeller  blade  outline  is  superimposed  on  Figure  4  with 
the  point  P  coinciding.  The  stations  between  the  line  piassing  P  and  the  leading  edge,  and 
between  P  and  the  trailing  edge  at  various  radial  stations,  including  the  stations  at  the  leading 
and  the  trailing  edges,  are  the  field  points  at  which  the  vortex  elements  are  assumed  to  be 
distributed. 

Then  the  computation  of  the  induced  velocity  components  at  point  P  (r^,  0^)  follows. 
This  is  done  by  computing  the  inte^ands  at  the  selected  chordwise  stations  and  integrating 
them  along  the  chord.  The  values  of  the  integrations,  which  represent  the  contribution  of  the 
vortex  systems  at  r  to  the  induced  velocities,  are  stored.  It  might  be  instructive  to  note  here 
that  in  the  calculation  of  Equations  [16c)  and  [16d)  the  strength  of  the  free  vortex  system, 
which  are  the  factors  enclosed  in  braces  in  these  equations,  may  be  represented  graphically 
as  shown  in -Figure  6.  They  may  be  considered  to  consist  of  two  parts,  one  is  the  contribution 
due  to  lifting  surface  and  the  other  due  to  the  superimposed  lifting  line. 

The  procedure  is  repeated  for  the  next  0^  as  point  P  moves  along  the  chord  toward  the 
trailing  edge  at  r^,  and  for  another  r,  and  so  on.  Then,  the  radial  integration  for  various  0^ 

•is  done  and  the  values  are  stored. 

When  the  calculations  for  the  tip-  and  hub-zones  of  the  “first  blade,”  and  for  the  other 
blades  are  completed,  the  calculation  for  the  mid-zone  of  the  first  blade  at  various  0q  follows. 

The  total  induced  velocity  components  at  P  (r^,  0^)  are  obtained  by  adding  the  values 
computed  for  the  mid-zone  to  the  radial  integrated  values  previously  computed.  The  entire 
procedure  is  then  repeated  for  the  next  and  so  on. 

The  downwash  0^)  at  each  point  P  specified  is  then  computed  from  the  total 

induced  velocity  components  according  to  Equation  [20].  The  downwash  at  points  where  mean 
line  ordinates  are  desired  is  extrapolated  or  interpolated  by  using  Equation  [21].  Finally,  the 

/  X 

ordinates  of  the  mean  line  at  any  point  I  r^,  - 

Basically,  the  new  computer  program  is  similar  to  the  original  program  in  the  method 
of  approach.  But  there  are  two  areas  in  which  improvements  have  been  made  in  this  new 
program  that  are  worth  mentioning.  One  is  in  the  method  of  setting  up  the  network  of  field 
points.  The  other  is  in  the  integration  of  the  induced  velocity  in  the  tip-  and  hub-zones. 

As  noted  above,  the  network  of  field  points  is  mapped  so  that  it  is  common  to  all  of 
the  control  points,  and  the  field  points  are  located  so  that  they  are  closely  spaced  in  the 


j  are  calculated  based  on  Equation  [22]. 
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neighborhood  of  the  control  point  and  are  symmetrical  about  it.  This  feature  of  symmetry 
improves  the  accuracy  of  numerical  computation,  especially  in  the  neighborhood  of  the 
singular  point. 

In  the  original  program  the  lifting  surface  is  divided  into  strips  by  cylinders  of  various 
radii.  Within  each  of  the  strips  the  pitch  distribution  and  the  circulation  functions  are  con¬ 
sidered  to  be  independent  of  the  radial  coordinate-  Thus,  the  integration  with  respect  to  the 
radial  coordinate  can  be  carried  out  first  by  using  the  mean  values  of  these  quantities  and 
then  the  chordwise  integration  is  done  numerically,  whereas  in  the  new  program  this  simplifi¬ 
cation  is  made  only  for  the  narrow  midstrip.  The  double  integrations  for  the  remaining  part 
of  the  lifting  surface  are  carried  out  numerically.  This  again  is  considered  to  be  an  improve¬ 
ment  for  better  accuracy  in  the  evaluation  of  the  integrals. 

RESULTS  OF  SAMPLE  CALCULATIONS 

The  computer  program  discussed  in  the  previous  section  has  been  completed.  A  number 
of  cases  have  been  worked  out  by  using  this  new  program  and  the  results  are  reported  herein. 

NUMBER  OF  BLADES 

The  first  case  is  a  5-bladed  propeller  similar  to  Propeller  3916,  reported  in  Reference  1, 
except  that  the  calculation  is  done  by  the  new  computer  program.  This  propeller  is  designated 
as  Propeller  3916A.  The  calculated  induced  velocities  and  the  camber  distributions  are  shown 
in  Table  1  and  in  Figure  15  (Appendix  C),  and  the  maximum  mean  line  ordinates  due  to  all  five 
blades  (Curve  1)  and  due  to  the  one  blade  (Curve  2)  are  plotted  in  Figure  7.  For  comparison, 
the  results  for  five  blades  reported  in  Reference  1,  calculated  from  the  original  computing 
program,  is  also  shown  (Curve  3).  The  difference  between  Curves  1  and  3  is  due  to  the  re¬ 
finement  made  in  the  present  computing  program. 

BLADE  OUTLINE 

It  is  worth  noting  that  the  camber  distribution  is  affected  by  the  shape  of  the  blade  out¬ 
line  selected.  When  the  blade  outline  of  Propeller  3916A  is  slightly  modified,  keeping  the 
same  expanded  area  ratio,  as  shown  in  the  following  table,  a  somewhat  different  camber 
distribution  is  obtained.  The  maximum  mean  line  ordinates  are  shown  also  in  Figure  7  (Curve  4). 
This  propeller  with  the  modified  blade  outline  is  designated  as  Propeller  3916B. 
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Original 

Modified 

0.2 

0.308 

0.294 

0.3 

0.372 

0.358 

0.4 

0.431 

0.430 

0.5 

0.462 

0.488 

0.6 

0.464 

0.488 

0.7 

0.435 

0.458 

0.8 

0.370 

0.364 

0.9 

0.268 

0.228 

1.0 

0.040 

0.000 

COMPARISON  WITH  TWO-DIMENSIONAL  CASE 

The  second  case  is  for  a  propeller  having  a  rectangular  expanded  outline  with  an  aspect 
ratio  of  10.  It  has  zero  hydrodynamic  pitch  ratio  and  a  uniform  bound-circulation  distribution 
along  the  radius.  This  fictitious  propeller  resembles  a  two-dimensional  wing,  although  the 
motion  involved  here  is  a  rotational  one  rather  than  a  linear  translational.  The  purpose  of  this 
calculation  is  to  determine  the  mean  lines  at  various  radii  induced  by  the  uniform  bound- 
circulation  distribution  only  and  to  compare  them  with  a  two-dimensional  airfoil.  The  calcu¬ 
lated  mean  line  distributions  at  radii  of  0.3,  0.6,  and  0.9  are  plotted  in  Figure  8  in  which  the 
NACA  fl  =  1  mean  line  distribution  is  also  shown. 

In  this  fictitious  case,  the  mean  lines  at  various  radii  are  quite  close  to  the  two- 
dimensional  case.  However,  for  the  actual  propeller,  the  means  lines  at  various  radii  obtained 
for  uniform  chordwise  loading  distribution  are  somewhat  different  from  the  c  =  1  mean  line  as 
shown  in  Figure  9  for  Propeller  3916.A.  Figure  10  shows  the  camber  distributions  of  Propeller 
3916A  at  various  radii  calculated  by  using  the  new  method  compared  with  those  calculated 
by  using  the  method  outlined  in  Reference  2  for  c  =  1  mean  line  with  camber  correction  factors 
Aj  and 


Camber  ratio,  corrected  =  0.05515  Aj  C ^ 
where  C,  is  lift  coefficient. 

L  f 

SKEWED  BLADES 

To  further  demonstrate  the  effect  of  blade  outline  on  camber  distribution,  a  third  case 
was  worked  out.  This  case  is  also  a  5-bladed  propeller  (Propeller  3917A).  The  loading, 
pitch,  and  blade  chord  length  distributions  are  the  same  as  those  for  Propeller  3916A,  except 
that  the  blades  are  skewed.  The  amount  of  the  skew  is  the  same  as  that  for  Propeller  3917.^ 
The  results  are  shown  in  Table  2,  and  the  comparative  camber  distributions  calculated  on  one 
blade  at  various  radii  are  plotted  in  Figure  11.  From  these  curves,  we  notice  that  blade  skew 
not  only  distorted  the  mean  line  but  also  changed  the  effective  angle  of  attack.  The  range  of 
angle  of  attack  for  this  particular  sample  problem  is  in  the  order  of  -0.4  degree  near  the  hub 
and  -0.6  degree  at  the  tip. 
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EXPANDED  AREA  RATIO 


To  show  the  dependence  of  camber  distribution  on  expanded  area  ratio  (EAR),  calcu¬ 
lations  were  made  for  propellers  having  the  same  basic  nlade  outline  as  Propeller  3916A  but 
with  different  EAR  ranging  f.X)m  0.5  to  1.2.  The  results  based  on  one  blade  are  shown  in 
Table  3  and  in  Figure  12.  In  this  figure,  the  maximum  mean  line  ordinates  due  to  the  bound 
circulation  only  are  also  shown.  The  contribution  due  to  the  other  four  blades  as  percentage 
of  the  contribution  due  to  one  blade  is  shown  in  Figure  13  for  the  various  EAR. 

LOAD  DISTRIBUTION 

Calculations  were  made  for  a  propeller  having  the  same  blade  outline  as  Propeller 
3916A  but  with  a  different  loading,  i.e.,  the  radial  bound  circulation  distribution  is  different. 

The  radial  load  distribution  is  altered  by  arbitrarily  changing  the  hydrodynamic  pitch  ratio  to 
a  constant.  This  constant  hydrodynamic  pitch  propeller  is  designates  as  Propeller  3916C. 

The  load  and  pitch  distributions  for  these  two  propellers  are  shown  in  Figure  14.  The  camber 
distributions  are  shown  in  Table  4.  The  comparative  camber  distributions  at  various  radii  are 
shown  in  Figure  15. 

DISCUSSION 

The  preceding  sample  calculations  were  presented  to  provide  a  clearer  picture  of  the 
information  that  can  be  obtained  from  the  program  and  of  the  various  parameters  that  affect 
the  mean  line  distribution. 

From  the  results  of  these  limited  sample  calculations,  we  may  make  the  following 
remarks : 

For  a  uniform  chordwise  load  distribution,  as  in  the  case  of  Propeller  39i6.A,  the  three- 
dimensional  mean  line  distribution  differs  somewhat  from  that  obtained  from  the  two-dimensional 
method  with  correction  factors. 

The  mean  line  distribution  of  blades  generally  depends  on  the  pitch  ratio,  radial  load 
distribution,  blade  chord  length  distribution,  and  blade  outline.  Consequently,  the  expanded 
area  ratio  and  the  skew  are  very  important  parameters  in  the  design  of  mean  line,  especially 
the  latter,  for  the  skew  not  only  distorts  the  mean  but  also  changes  the  effective  angle  of 
attack. 

The  presence  of  other  blades  has  some  effect  upon  the  mean  line  of  a  blade.  It  has 
more  effect  in  the  hub  region  than  it  has  near  the  tip.  The  contribution  to  the  mean  line  due 
to  other  blades  also  increases  as  the  expanded  area  ratio  increases  for  propellers  having  the 
same  blade  outline. 

It  may  be  of  interest  to  mention  that  some  of  the  these  findings  disagree  with  the  con¬ 
clusions  drawn  in  Reference  3  in  which  van  Manen  and  Bakker  reported  the  numerical  results 
of  a  number  of  propeller  design  problems  using  lifting-surface  theory.  They  reported,  in  part. 
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that  for  a  uniform  pressure  distribution  along  the  chord,  the  three-dimensional  camber  line  of 
a  propeller  blade  section  is  identical  to  the  camber  line  of  the  two-dimensional  profiles  and 
that  the  effect  of  radial  load  distribution  is  not  very  pronounced.  They  also  reported  that  the 
effect  of  skew  back  is  negligible.  Based  on  the  evidence  shown  in  this  report,  their  conclu¬ 
sions  appear  to  be  doubtful. 


FUTURE  PLANS 

To  obtain  some  experimental  verification  of  the  theory,  tests  of  model  propellers  de¬ 
signed  by  using  this  method  are  contemplated.  The  tests  are  to  be  done  in  conjunction  with 
other  research  projects  that  are  currently  going  on  at  the  David  Taylor  Model  Basin.  It  is 
also  felt  that  the  test  of  the  theory  could  be  easily  accomplished  by  comparing  the  results 
from  the  inverse  problem  of  performance  prediction  for  the  steady  inflow  case  with  reliable 
published  open-water  experimental  data.  If  close  agreement  between  the  theoretical  and 
experimental  results  is  obtained,  of  course,  the  validity  of  the  theory  is  firmly  established. 
If,  however,  some  discrepancies  are  found,  we  hope  to  be  able  to  derive  some  consistent 
parametric  relationships.  When  this  is  accomplished,  propeller  open-water  design  charts 
can  be  constructed  with  data  obtained  from  computations  for  a  large  range  of  parameters 
rather  than  from  costly  experiments. 

With  these  considerations,  we  shall  proceed  with  the  work  on  the  remaining  part  of  the 
overall  program.  The  work  on  the  case  of  an  arbitrary  chordwise  load  distribution  is  already 
underway  and  will  be  reported  in  the  near  future.  It  is  anticipated  that  the  work  on  the  in¬ 
verse  problem  of  predicting  performance  for  a  given  propeller  design  will  commence  shortly. 

Also  anticipated  is  the  inclusion  of  the  thickness  effect  in  the  investigation.  Work 
in  this  area  has  been  done  recently  by  Kerwin  and  Leopold;  their  results  are  reported  in 
References  4  and  5.  Their  work  may  be  incorporated  in  the  present  computation  to  take 
account  of  the  thickness  effect. 

Finally,  from  the  investigation  of  the  unsteady  inflow  case,  it  is  hoped  that  one  can 
estimate  the  time  dependent  force  fluctuations  produced  by  a  propeller  operating  in  a  nonuni¬ 
form  wake  field  behind  a* ship.  Thus,  the  phenomena  of  cavitation  due  to  unsteadiness  and 
propeller-induced  vibration  and  noise  could  be  better  understood. 
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Figure  2  —  Coordinate  System  in  Projected  Plane  and  Regions  of  Integration 
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Figure  3  —  \'elocity  Diagram  at  Radius 


Figure  5  —  Field  Points  for  Point  P  on  Blade 
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Figure  6  -  Graphic  Representation  of  Free  Vortex  Strength  Distribution  Along  a  Helical  Line  at  r  Due  to 

a  Lifting  Surface  and  a  Superimposed  Lifting  Line  at  <;S  =  0 


MAXIMUM  MEAN  LINE  ORDINATES,  (Cm)  IN  INCHES 
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Figure  7  -  Maximum  Mean  Line  Ordinates  for  Propellers  3916,  3916A,  and  3916B 


Figure  8  -  Camber  Line  Distribution  of  a  Two-Dimensional  Propeller  Blade 
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PERCENT  or  MAXIMUM  CAMOSR 
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Figure  9  —  Camber  Line  Distribution  of  Propeller  3916A 
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MEAN  LINE  ORDINATES,  C^,  IN  INCHES 
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DISTANCE  FROM  LEADING  EDGE  AS  FRACTION  OF  CHORD,  x/L^ 

Figure  10  -  Comparison  of  Camber  Distribution  for  Propeller  3916A 


28 


DISTANCE  FROM  LEADING  EDGE  AS  FRACTION  OF  CHORD,  x/Uj. 


Figure  11  —  Comparison  of  Camber  Distributions,  Symmetrical  versus  Skewed  Blade  Propellers 
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MEAN  LINE  ORDINATES,  (Cm)n^ax  In  Inches 
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Figure  12  —  Maximum  Mean  Line  Ordinates  for  Propellers  with  Different  EAR 
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Figure  13  —  Percentage  Contribution  to  Ma^cimum  Moan  Line  Ordinates  Due 
to  Other  Blades  for  Propellers  with  Different  EAR 
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Figure  14  —  Pitch  Ratio  and  Bound  Circulation  Distribution  for  Propellers  3916  and  3916C 
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MEAN  LINE  ORDINATES,  C,n,  IN  INCHES 


0.10 


0.09 


-  PROPELLER  39 16A 

- PROPELLER  3916C  (with  Constant  Hydrodynamic  Pitch) 

Calculation  based  on  5  Blades 


DISTANCE  FROM  LEADING  EDGE  AS  FRACTION  OF  CHORD,  x/Lj 

Figure  15  —  Comparison  of  Camber  Distribution  for  Propellers  3916A  and  3916C 
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TABLE  1 

Calculated  Induced  Velocities  and  Camber  Distribution  for  Propeller  3916A 
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0*0035 

0*0056 

0*0069 

0*0077 

0*0080 

0*0004 

0*0005 

0*0009 

0*0017 

0*0029 

0*0047 

0*0059 

0*0066 

0*0069 

NOTES:  1.  CalculMloa  bMod  on  one  bfjde  ooly. 
3.  CftlculAted  cAmbor  in  tnebeo. 


o.osoo 

0*)  164 

0.0036 

0.0736 

0.0S6I 

0.0360 

0*0260 

O.Otll 

0*0026 


O.OS60 

•0.0343 

•0.027S 

•0*0217 

•0*0170 

•0*0111 

•0*0073 

•0*0033 

•0.0606 


0*9900 

0.1213 

0*0979 

0*0769 

0*0609 

0*0396 

0*0360 

0*0116 

0*0027 

0*9900 

0*6976 

0.745*' 

1*4943 


0*0005 

0*0006 

0.0012 

0*0022 

0*0040 

0*0068 

0*0066 

0*0100 

0.0105 


0*0003 

0*0004 

0*0006 

0*0015 

0*0027 

0*0046 

0*0059 

0*0067 

0*0070 
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TABLE  2 


Calculated  Induced  Velocities  and  Camber  Distribution  for  Propeller  3917A  (Skewed) 

IX^UT  C>TA 

It 

0*00390  0*IOC00  2*OOOOOt90*OOCOO  / 

0*00300  0*01000  0*01600  0*02«00  0*03300  0*03000  0*07000  0*00600  0*12000 

0*17000  0*22000  0*26000  0*33300  0*44C00  0*33200  0*70000  0*04300  t*<l0000 

3  0  0  0  0  -0  -0 

0*01230  0*10166  0*23724  0*40626  0*30172  0*76276  0*69634  0*08014 

0*01230  0*02300  0*03000  0*10000  0*20000  0*40000  0*60000  0*60000  0*93000 

0*20000  O«3C00O  0*40000  0*30000  0*60000  0*70000  0*60000  0*90000  1*00000 

0*12300  0*20900  0*26400  0*26600  0*27300  0*22700  0*14900  0*03300-0*11300 

0*30600  0*37200  0*43000  0*46200  0*46400  0*43300  0*37000  0*26600  0*04000 

1*66210  1*32633  0*90766  0*70429  0*36622  0*46392  0*36696  0*32722  0*27960 


0*  0*01730  0*02266  0*02426  0*02263  0*01936  0*01463  0*00696-0* 

0*23000  0*30000  0*40000  0*30000  0*60000  0*70000  0*60000  0*90000  0*93000 

1*00000  1*00000  1*00000  1*00000  1*00000  1*00000  1*00000  1*00000  1*00000 

6*14300  0*16902  0*20661  0*22173  0*20763  0*17661  0*13364  0*08937  0*06300  UftCl’o) 

6*16010  0*16676  0*16039  0*13704  0*10317  0*07463  0*04669  0*02767  0*01770 

1*00600  0*62600  6*00000  1*00000  1*00000  Card  9 


Card  1  B  for  explanatlco) 

Card  2 
Card  3 

Card  5 
Card  f> 

Card  7 

Card  B  r 

L/^r)/D 
Lt(r>/D 
tao  3i(7) 


CUT6CT 

6XUt  l»»0UCf0  VCLCCITV  C0P6C6eMT 


6X/N0 

0*2300 

O*3C00 

0*4000 

0*3000 

0*6000 

0*7000 

0*6000 

0*9000 

0*9500 

X/LT 

0*0123 

0*0609 

0*1113 

0*1136 

0*1427 

0*1442 

6*1770 

0*1416 

0*1922 

0*1460 

0*1017 

0*0373 

0*0669 

0*0604 

0*0713 

0*0660 

0*0472 

0*0964 

0*1237 

0*1343 

0*2372 

0*0430 

0*0400 

0*0232 

0*0200 

0*0320 

0*0430 

0*0460 

0*0644 

0*1634 

0*4063 

0*0341 

0*0163 

-0*0067 

-0*0122 

-0*0049 

0*0090 

0*0213 

0*0336 

0*0720 

0*3917 

0*0226 

-0*0063 

-0*0322 

-0*0423 

-0*0333 

-0*0201 

-0*0067 

0*0233 

0*0333 

0*7628 

0*0126 

-0*0310 

-0*0373 

-0*0716 

-0*0630 

-0*0460 

-O-tOJSO 

-0*0036 

-6*0003 

0*6963 

-0*0133 

-0*0369 

•0*0676 

-0«I049 

-0*0973 

-0*0763 

-0*0420 

-0*0333 

-0*0323 

0*9601 

-0*0429 

-0*1002 

-0*1296 

-0*1312 

-0*1446 

-0*1231 

-0*1031 

-0*0729 

-0*0927 

tAN«C6TIAL  INCUCCC  VCUOCITT 

CCRPONCNT 

6X/60 

0*2500 

0*3600 

0*4000 

0*3000 

0*6000 

0*7000 

0*6000 

0*9000 

0*9906 

2/UT 

0*0123 

-0*1953 

-0*2042 

-0*1437 

-0*1423 

-0*1224 

-0*6994 

-0*0756 

-0*0462 

-0*6479 

0*1017 

-0*1479 

-0*1416 

-0*1073 

-0*0690 

-0*0716 

-0*0347 

-0*0416 

-0*0436 

-0*6373 

0*2372 

-0*1066 

-0*0934 

-0*0632 

-0*0460 

-0*0374 

-0*0266 

-0*0201 

-0*0269 

-6*6279 

0*4063 

-0*0362 

-0*0372 

-0*6133 

-0*0069 

-0*0032 

-0*0043 

-0*0064 

-0*0179 

-0*0179 

0*3917 

-C«0022 

0*0231 

0*0346 

0*0333 

0*0243 

0*0139 

0*0073 

-0*9064 

-0*6646 

0*7626 

0*0463 

0*0794 

0*0606 

0*0703 

0*0311 

0*0323 

0*0206 

0*0046 

0*6648 

6*6983 

0*0952 

0*1291 

6*1217 

««1036 

6*0743 

0*0864 

0*0334 

6*6193 

6*6144 

6*9801 

0*1404 

0*1793 

6*1683 

0*1404 

0*1682 

0*0732 

0*0801 

0*0269 

6*6326 

titouceo  ooBNuAs^*  cue  rc  ai.l 

6CA066 

■x/no 

0*2300 

0*3000 

0*4006 

9*3000 

0*4000 

0*7006 

0*8600 

6*9609 

6*8366 

A/LT 

6*0123 

0*2073 

0*2266 

0*1936 

0*1966 

0*2033 

0*2024 

0*1969 

0*2C33 

6*1747 

0*1017 

0*1946 

0*1822 

0*1166 

0*1093 

0*1030 

0*1621 

0*1049 

6*1330 

6*1398 

0*2372 

0*1133 

0*0977 

0*0397 

0*0440 

0*0442 

0*0311 

0*0320 

0*0692 

0*1072 

0*4063 

0*0657 

0*0^91 

0*0033 

-0*0060 

-0*0017 

0*0106 

0*0221 

0*0963 

0*074t 

0*3917 

0*0143 

-0*02^20 

-0«047i 

-0*0340 

-0*0429 

-0*0241 

>0*0089 

0*0260 

0*0336 

0*7626 

-0*0332 

-0*0612 

-0*0966 

-0*0992 

-0*0816 

-0*0372 

-0*0390 

-0*0030 

-0*0017 

0*6963 

-0*06eS 

-0*1373 

-0*1466 

-0*1433 

-0*1224 

-0*0933 

-0*009^ 

-0*0364 

-6*0351 

0*9601 

-0*1403 

-0*2023 

-0*2071 

-0*2044 

-0*1777 

-0*1423 

-0«||43 

-0*0763 

-0*0961 

HX/NC 

0*2300 

0*3000 

0*4000 

0*3000 

0*6000 

0*7000 

0*6000 

0*9000 

0*9806 

P/0 

1*1360 

1*1377 

1*1406 

1*IC63 

1*0673 

1*0202 

0.9726 

0*9252 

0*6676 

L*b* 

1*3624 

1*6720 

2*1120 

2,2660 

2*1040 

1*6140 

1*1920 

0*2640 

-0*3133 

CHCNO 

2*7213 

2*9760 

3*4400 

3*6960 

3*7120 

3*^t00 

2*9600 

2*1440 

1*4943 

CA6M6  OISTRteOTICK 

X/LT 

0*0123 

0*0032 

0*0036 

0*0046 

0*0043 

0*0040 

0*0033 

0*0024 

0*0016 

0*0009 

0*0230 

0*010  1 

0*0112 

0*0093 

0*0086 

0*0674 

0*0042 

0*0046 

0*0031 

0*0017 

0*0300 

0*0194 

C*0210 

0*0171 

0*0136 

0*0136 

0*0112 

0*0064 

0*0097 

0*0034 

0*1000 

0*0337 

0*0373 

0*0299 

0*0270 

0*0232 

0*0187 

0*0143 

0*0102 

0*0043 

0*2000 

0*0626 

0*0630 

0*0476 

0*0411 

0*0331 

0*0264 

0*0210 

0*0170 

0*0114 

0*4000 

0*1004 

0*0926 

0*0624 

0*0497 

0*0433 

0*0372 

0*0293 

0*0264 

0*0169 

0*6000 

0*1160 

0*0962 

0*0344 

0*039A 

0*0367 

0*0333 

0*0306 

0*0319 

0*0233 

0*6000 

0*1100 

0*0720 

0*0248 

0*0109 

0*0144 

0*0243 

0*0234 

0*0323 

0*0241 

0*9300 

0*0673 

0*0354 

-0*0144 

-0*0246 

-0*0093 

0*0064 

0*0164 

0*0294 

0*0223 

CAPeCR  RAtIO 

X/tT 

0*0125 

0*0019 

0*0020 

0*0014 

0*0012 

0*0011 

0*0009 

0*0006 

0*0007 

0*0006 

0*0230 

0*0037 

0*0038 

0*0027 

0*0023 

0*0020 

0*0014 

0*0013 

6*0014 

0*0012 

0*0300 

0*0071 

0*0071 

0*0030 

0*0043 

0»0037 

0*0032 

0*002S 

0*0027 

0*0023 

0*1000 

0*0131 

0*0126 

0*0067 

0*0073 

0*0043 

0*0034 

0*0046 

0*0046 

0*0042 

0*2000 

0*0231 

0*0212 

0*0136 

0*0111 

0*0093 

0*0062 

0*0074 

0*0079 

0*0074 

0*4000 

0*0369 

0*0312 

0*0161 

0*0134 

0*0917 

0*0107 

0*0099 

0*0123 

0*0127 

0*6000 

0*0426 

0*0323 

0*0138 

0*0107 

0*0099 

0*0102 

0*0103 

0*0147 

0*0184 

0*6000 

0*0404 

0*0243 

0*0072 

0*0030 

0*0043 

0*0071 

0*0066 

0*0190 

0*0142 

0*9300 

0*0321 

0*0119 

-0*0042 

-0*0067 

-0*0023 

0*0023 

0*0048 

0*0137 

4*4146 

NOTES:  t.  Calculation  bated  on  five  bUdet. 
2.  Calculated  camber  In  Inches. 
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TABLE  3 

Camber  Distribution  for  Propellers  with  Different  Expanded  Area  Ratios 


EAR  •  0. 50 


»)X/RO 

O*2S00 

0*3000 

0*4000 

3*5000 

8*6000 

0*7000 

0*8000 

0#9000 

0*9500 

P/O 

]*t500 

1*1577 

1*1406 

1  *  1063 

1*0673 

1 *0202 

0*9726 

0*9252 

0*8976 

0*7tt9 

0*7827 

0*9068 

0*9720 

0*9763 

0*9152 

0*7785 

0*5639 

0*3920 

Chord 

1*43)4 

1*5654 

1*80  >4 

1 *9441 

1*9525 

1 *8305 

1*5570 

1* 1277 

0*7860 

CAMtlCP  OISTRtOUTION 

X/LT 

0*0009 

0*0022 

0*0020 

0*0028 

9*00?6 

0*0021 

0*0015 

0*0010 

0*0006 

0.0003 

0*0I2S 

0*0027 

0*0035 

0*0035 

0*0033 

0*0026 

0*00)9 

0*0013 

0*0007 

0*0004 

0*0200 

o*ooso 

0*0065 

0*0065 

0*0061 

0*0049 

0*0036 

0*0024 

0*001 3 

0*0007 

0*0S00 

0*007t 

0*01  16 

0*0116 

0  *  0  1  0;8 

0*0088 

0*0065 

0*0043 

0*0024 

0*0014 

0*1000 

3*0157 

0*0108 

0*01 >5 

0*0181 

0*0148 

0*0109 

0*0074 

C*0042 

0*0026 

0*2000 

0*0257 

0*0321 

0*0309 

0*0288 

0*02  15 

0*0172 

0*011  7 

0*0067 

0*0043 

0*3000 

0*0325 

0*0402 

0*0302 

0*0357 

0*0289 

0*0213 

0*0146 

0*0084 

0*0054 

0*4000 

0«0365 

0*0449 

0*04<^5 

0*0397 

0*0322 

0*02  36 

0*0162 

0*0094 

0*0061 

0*9000 

0*0379 

0*0463 

0*0438 

0*0410 

0*0332 

0*0244 

0*0167 

0*0097 

0*0064 

EAR  -  0.75 


RX/RO 

0*2500 

0*3000 

0*4000 

0*5000 

0*6000 

0*7000 

0*8000 

0*9000 

0*9500 

P/0 

1  *  1560 

1* 1577 

1*1405 

I* 1063 

1*0673 

1*0202 

0.«/726 

0*9252 

0*8976 

1*0664 

1 *1725 

1*3585 

1*4562 

1*4625 

1*371  t 

1 • 1662 

0*8447 

0*5872 

CHORD 

2*1444 

2*3451 

2*7107 

2*9124 

2*9251 

<*7422 

2.3)25 

1 *6875 

1*1775 

CAMSCR  OISTRieuriON 

x/u? 

0*0099 

0*0025 

0*0032 

0*0032 

0*0030 

0*0025 

0*0018 

0*0013 

0*0007 

0.0004 

0*0125 

0*0030 

0*0019 

0*0039 

0*00)7 

0*0030 

0*0022 

0*0016 

0*0009 

0*0005 

0*0250 

0*0058 

0*0074 

0*0074 

0*0069 

0*0057 

0*0042 

0*0029 

0*001  7 

0*0010 

0*0500 

0*0106 

0«0I55 

0*01 32 

0*0124 

0*0102 

0*0076 

0*0052 

0*0031 

0*0019 

0*1000 

0*0185 

0«02 35 

0*0224 

0*0212 

0*0174 

0*0128 

0*0088 

U0053 

0*0033 

0*2000 

0*0307 

0*0385 

0*03M 

0*034  1 

0*0277 

0*0200 

0*0l4l 

6*0085 

0*0057 

0*3000 

0*0391 

0*0405 

0*0450 

0*0424 

0*0344 

0*0251 

0*0176 

U.0107 

0*0073 

0*4000 

0*0442 

0*0546 

0*0501 

0*0472 

0.03U3 

0*0279 

0*0195 

0*0120 

0*0083 

0*5000 

0*0460 

0*0565 

0*0516 

0*0409 

0.0395 

0.0288 

0*0202 

0*0124 

0*0066 

EAR  >  1.20 


RX/RO 

0*2900 

0*3000 

0*4000 

0*9000 

0*6000 

3*7000 

0*8000 

0*9000 

0*9500 

P/0 

1*1560 

1*1577 

1*1406 

1*1063 

1*0673 

1*0202 

0*9726 

0*9292 

0*8976 

1*7052 

1*8749 

2*1722 

2*3265 

2*3386 

2*1924 

1  *6648 

1 *3507 

0*9389 

CHORD 

3*4289 

3*7498 

4*3344 

4.6570 

4*6771 

4*3848 

3*  7296 

2*7014 

1*0829 

CAMtCR  OISTRIDUTION 

X/LT 

0*0099 

0*0030 

0*0040 

0*0038 

0*0038 

0*0032 

0*0024 

0*0017 

0*0010 

0*0006 

0*0129 

0*0038 

0*0049 

0*00  > 

0*0047 

0*0039 

0*0030 

0*6022 

0*0013 

0*0007 

0*0290 

0*0072 

0*0094 

0*0090 

0*0086 

0*0073 

0*0057 

0*0040 

0*0024 

0*0014 

0*090C 

0*0133 

0*0)72 

0*0163 

0*0157 

0*0128 

0*0102 

0*0072 

0.00** 

0*0027 

0*1000 

0*0238 

0*0309 

0*0283 

0*0271 

0*0219 

0*0170 

0*0122 

0*0075 

0*0050 

0*2000 

0*0403 

0*0912 

0*0482 

0*0442 

0*0359 

0*0769 

0*0194 

0*0122 

0*0005 

0*3000 

0*091 7 

0*0693 

0*0979 

0*0555 

0*0450 

0.0337 

0*0241 

0*0153 

0*0110 

0*4080 

0*0986 

0*0737 

0*0647 

0*0621 

0*0903 

0*0379 

0*0267 

C',01  72 

0*0125 

0*90C«0 

0*0613 

0*0764 

0*0687 

0*0643 

0*0920 

0*0388 

0*0277 

0*0179 

0*0131 

NOTES:  1.  Blade  outline  of  Propeller  3916A  which  haa  an  EAR  •  0.95  la  uaed  aa  baalc  blade  otiUlne. 

2.  Calculatlona  are  baaed  on  one  blade  only. 

3.  Calculated  camber  distributions  are  In  Inches. 
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TABLE  4 


Calculated  Induced  Velocities  and  Camber  Distribution  for  Propeller  3916C 


tRRVT  CATA 

IS 

C»t6  1 

(ttr  Appendix  B  for  cxpUnxtlon) 

••00900 

i  0*10000  2*00000190*00000 

C«rd  2 

••oosoo 

0*01000 

0*01400 

0*02400 

0*09960 

0*09000 

0*07000 

0*09600 

0*12900 

C4rd  3 

••17000 

0*22000 

0*2«000 

0*99200 

0*44000 

0*99200 

0*70900 

0*94600 

1*00000 

% 

• 

9 

1 

9 

-0 

-0 

Card  } 

e*ooff j 

0*09099 

0*lia62 

0*20414 

0*29966 

0*96196 

0*44917 

0.49C07 

Card  6 

••OOffJ 

0*01290 

0*02900 

0*09000 

0*10000 

0*20000 

0*90000 

0*40000 

0*90000 

Card  7 

0*20000 

0*90000 

0*40000 

0*90000 

0*60000 

0*70000 

0*60000 

0*90000 

1*00000 

Card  8 

r 

0*19400 

0*l««00 

0.21990 

0*29100 

0*29200 

0*21790 

0*16900 

0*13400 

0*02000 

l.l(r)/D 

0*90000 

0*37200 

0*43000 

0*46200 

0*46400 

0*49900 

0*97000 

0*26600 

0*04000 

lt(r)/B 

l*97US 

1*04777 

0*7e9a9 

0*62666 

0*92966 

0*4*904 

0*99291 

0*94926 

0*91499 

can  ®j(r) 

0* 

0*01067 

0*01491 

0*01696 

0*016*1 

0*01697 

0*01699 

0*01997 

0* 

r(r> 

0*29000 

0*90000 

0*40000 

0*90000 

0*60000 

0*70000 

0*60000 

0*90000 

0*99000 

1*00000 

1*00000 

1*00000 

1*00000 

1.00000 

1*00000 

1*00000 

1.00000 

1.00000 

Vx(ro) 

0*07900 

0*00299 

0*12094 

0*19997 

0*19267 

0*16219 

0*16690 

0*17973 

0*17900 

“.(Co) 

0*00990 

0*00799 

0*09910 

0*06617 

0*06041 

0*07306 

0*06661 

0*06902 

0*06000 

“t<Co> 

uooooo 

0*a2600 

•*00000 

1  *00000 

1*00000 

Card  9 

CUT^T 

AXIAL  INOVICCO  VCLOCfTT  COMPONCNT 

RX/RO  0*2900  0*9000 

0*4000 

0*9000 

0*6000 

0*7000 

0*6000 

0*9000 

0*9900 

X/LT 

0*0099 

0*021 4 

0*0999 

0*0649 

0*0861 

0*1069 

0*1946 

0*1777 

0*2487 

0*2667 

0*0906 

0*0104 

0*0970 

0*0969 

0*0996 

0*0667 

0*0669 

0*1209 

0*1602 

0*2219 

0.1I46 

0,0090 

0*0246 

0*0292 

0*0946 

0*0410 

0*0979 

0*0606 

0*1967 

0*1732 

0*2041 

0*0007 

0*0199 

0*0126 

0*0209 

0*0267 

0*0961 

0*0999 

0*0993 

0*1970 

0*2999 

0*0009 

0*0071 

0*0064 

0*0109 

0*0111 

0*0229 

0*0392 

0*0672 

0*0991 

0*9614 

0*0002 

0*0092 

0*0022 

0*0096 

0*0081 

0*0121 

0*0199 

0*0987 

0*0964 

0*4492 

0*0007 

0*0022 

0*0007 

0*0024 

0*0099 

0*0091 

0*0069 

0*0164 

0*0296 

0*4901 

0*0010 

0*0009 

•O.OOOi 

0*0006 

0*0006 

0*0010 

0*0017 

0*0027 

0*0097 

tancchtial 

RX/RO 

INOUCeO  VELOCITT 

0*2500 

C0MP0N6NT 

0*9000 

0*4000 

0*9000 

0*6000 

0*7000 

0*6000 

0*9000 

0*9900 

X/LT 

0*0099 

•0*1014 

•0*1112 

•0*0996 

•0*0926 

•0*0699 

•0*0777 

•0*0766 

•0*0914 

•0*0900 

0*0906 

•0*0626 

•0*0699 

•0*0772 

•0*0700 

•0*0619 

-0*0996 

•0*0996 

-0*0670 

•0*0749 

0*1166 

•0*0674 

•0*0716 

•0*0602 

•0*0999 

•0*0441 

•0*0994 

•0.0989 

•0*0912 

•0*0987 

9*2041 

•0*091 7 

•0*0949 

•0*04*6 

•0*0964 

-0*0914 

•0*0269 

•0*0271 

-0.0973 

•0*0469 

0*2999 

•0*0999 

•0*0974 

•0*0901 

•0*0299 

•0*0194 

-0*0172 

•0*0174 

•0*0292 

-0*0919 

0.94I4 

•0*021 9 

•0*0220 

•0*0172 

•0*0144 

•0*0112 

•0*0099 

•0.0099 

•0*0l49 

-0*0197 

0*4492 

•0*0096 

•0*0099 

•0*0072 

•0*0062 

-0.0047 

-0*0040 

•0*0042 

•0*0061 

•O*0087 

0*490* 

•0*0029 

•0*0019 

-0*0012 

-0.0019 

•0*0009 

•0*0006 

-0.0009 

•0*0010 

•0*0020 

fNOilCCO  OaVMVASH  OUC  TO  ALL  OLAOCS 


RX/RO 

0*2500 

0*9090 

0*4000 

0*9000 

0*6000 

0.7000 

0*6000 

0*9000 

0*9900 

x/lt 

0*0099 

0.0927 

0*1190 

6* 1127 

0*1222 

0.1994 

0.1946 

0*1941 

0*2649 

0*2014 

0*0906 

0*0717 

a*09qi 

0*0709 

0*0644 

0*0693 

0.1036 

0*13/9 

0*19/2 

0*2336 

e.ilio 

0*0999 

0*6669 

0*0999 

0*0976 

0*0967 

0*0664 

0*0699 

0.1499 

0.1629 

0*2041 

9*0409 

0*0466 

0.0976 

0*0376 

0*0302 

0*0439 

0*0614 

0*1060 

0*1 446 

0*2999 

0*0269 

0*0920 

0.02J6 

0*0227 

0.0109 

0*0279 

0*0372 

0*0718 

0*0963 

0*3014 

0*0166 

0*0199 

0.0129 

0*CI29 

0.0123 

0*0190 

0*0222 

0*0413 

0.061  7 

0*4492 

0*0079 

0*0084 

0.0090 

0*0063 

0.0091 

0*0063 

0*0099 

0*0176 

0.0272 

0*4901 

0*0026 

0*0017 

0.6009 

0*0012 

0.0010 

0*0012 

0*0019 

0.00/7 

0.0060 

RX/RO 

0*2900 

0*9000 

0.4000 

0*5000 

0.6000 

0*7000 

0*6000 

0.9000 

0.9900 

P/0 

0*9679 

0*9679 

0.9879 

0*9079 

0.9879 

0*9679 

0*9679 

0*9679 

0.9879 

L*e* 

1*9999 

1*4660 

1.7240 

1*0460 

1.6960 

1*7400 

1.4600 

1*07/0 

0.7492 

CHORD 

2*721 9 

2*9760 

9*4400 

3*6960 

3.7120 

3*4000 

2*9600 

?• 14a0 

1.4943 

CAMIEA  OISTAteU) lOH 
X/LT 


0*0099 

0*0019 

0*0024 

0*0022 

0*0022 

0.0021 

0.0020 

0*0019 

0*001  7 

0.0011 

0*0129 

0*0024 

0*0090 

0*0026 

0*0026 

0*0026 

0.0029 

0*0023 

0*0021 

0.0014 

0*0290 

0*0049 

0*0096 

0*0092 

0*0062 

0*0049 

0.0046 

0*0044 

0*0039 

0.0028 

0*0900 

0*0063 

0*0106 

0*0099 

0.0092 

0*0067 

0.0062 

0*0076 

0*0071 

0.0063 

0*1000 

0*0147 

0*0186 

0.0166 

0.0199 

0*0143 

0*0136 

0*0131 

0*0122 

0.0096 

0*2000 

0*0247 

0*0906 

0*0261 

0.0241 

0*0217 

0*0200 

0*0204 

0*0201 

O.0162 

0*3000 

0*031 9 

6*0969 

0*0306 

0.0292 

0.0299 

0*0261 

0*0260 

0*0266 

0.021  I 

0*4000 

0*0397 

0.0432 

0*0939 

0.0920 

0.0260 

0*0279 

0*0276 

0*0286 

O.0241 

0*9000 

0*0972 

0*0449 

0*0346 

0.0929 

0.0266 

0*0203 

0.0266 

0*0297 

0.0261 

CAXeEX  RATIO 
X/LT 


0*0099 

0*0007 

0*0008 

0*0007 

0*0006 

0*0006 

0.0006 

0*0006 

0*0008 

0*0008 

0*0129 

0.0009 

0*0010 

0*0006 

0*0007 

0*0007 

0.0007 

0*0006 

0*0010 

0*0010 

0*0290 

O.OOt  7 

0*0019 

0*0019 

0*0014 

0*001  3 

0.0019 

0.0019 

0*0016 

0*0019 

0*0500 

0*0031 

0*0099 

0*0027 

0*0029 

0*0024 

0*0024 

O.C026 

0*0033 

0*0039 

O.I^CrO 

0*0094 

0*0062 

0*0046 

0*0042 

0*0096 

0.  10)9 

0.0044 

0*0097 

0*0064 

5*2000 

0*0091 

0*0103 

0*0079 

0*0069 

0*0099 

0*0060 

0*0069 

0*0094 

0*0109 

0*9000 

0*0116 

0*0129 

0*0090 

0.0079 

0*0070 

0*0072 

0*0064 

0*0119 

0*0141 

0*4000 

0*0191 

0*0146 

0*0099 

0.0067 

0*0079 

0*0079 

0.0093 

0*0134 

0*0161 

0*9000 

0*0197 

0*0161 

C.OIOI 

0*0069 

0.0078 

0*0061 

0*0096 

0*0136 

0.0166 

NOTES:  1  Calculation  baaed  on  liv«  bladea. 
2.  Calculated  camber  In  Inches. 
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APPENDIX  A 


EVALUATION  OF  INTEGRALS 


It  has  been  stated  in  the  text  that  the  integrals  of  Equation  [16]  have  a  singularity  at 
f  =  Tq  and  0  =  (or  ^  =  0).  The  integrals  are  integrated  numerically  with  this  singular  point 
isolated.  The  justification  for  this  procedure  is  that  the  contribution  to  the  final  mean  line 
due  to  the  small  region  that  includes  thq  singular  point  is  negligible.  To  show  that  this  is 
indeed  the  case,  we  shall  first  take  a  look  at  the  induced  velocities  and  then  their  contri¬ 
butions  to  the  final  mean  line. 


INDUCED  VELOCITIES 

The  integral  to  be  investigated  which  represents  the  contribution  of  the  small  area  to 
the  induced  velocity  components  is  in  the  following  general  form: 


where  Sq«1»  <q«1>  represents  the  loading  function,  and  0p,  r,  0)  is  the 

kernel  function  which  has  the  form 


V  ^)  = 


(r,  0) 


+  f2  _  2  Tjj  r  cos  (0  -  0^)  +  {p/nf  (0  - 

lim 


3/2 


To  facilitate  the  evaluation,  we  introduce  a  limiting  process  of  such  that 


lim 

^'’o’  X-*0 


.'b‘^^0  *0 


"o"^o  W^  +  ^o +r^-2ror  cos  (0-0p)  +  (pA)^(0-0jj)^] 


3/2 


Physically,  the  introduction  of  this  limiting  process  means  that  mathematical  evalua¬ 
tions  are  performed  over  an  area  on  a  surface  which  is  X  distance  away  from  the  control  point 
on  the  lifting  surface,  and  that  in  the  limit  the  fictitious  area  is  brought  infinitesimally  close 
to  the  control  point.  This  means  that  the  velocity  induced  on  the  surface  /Y  =  0  is  a  continuous 
function  of  X  including  X  =  0. 

With  this  procedure  in  mind,  we  shall  attempt  to  deal  with  each  of  the  four  integrals 
involved. 
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Equation  [16a] 


/o  +  ^o  r^o'^'o 


lim  'b  r  °  r''°  ■  y(')  -  ^o) 

4^  1  J - 

'b“^0  ^o“<0  iX^+r^  +  r^-%r  cos{e-d^)  +  {p/nf{e~dQf] 


For  a  small  <q,  we  may  substitute  O-Oq  for  sin  (O-d^)  and  f  -  ~  for 

cos  (0-0q),  and  for  a  linear  approximation  we  can  express 


y(r)  =  y  (<■())  +  yVo) 

t  =  0-6^,  l«l  «1 

5  =  r  -  fg,  |S1  «  1 


we  have 


^o)“  -Y-oi"  77 


r  ^  f*”  f*”  lHr„)H-y-(g81 

1  *’  -«„  -*0  lX^*SUa^t^r 


d*  dS 


where 


The  integration  of  the  second  term  with  respect  to  5  is  zero,  which  justifies  the 
omission  of  the  second  term  in  the  calculation  for  the  mid-zone,  and  the  integration  of  the 
first  term  also  vanishes  for  0^  -  €q>  0q>  0^  +  e^.  Therefore,  we  only  have  to  consider  the 
cases  of  0q>  -  <o  *  ‘o’  points.  Integrating  the  first  integral 

with  respect  to  <  and  assuming  0^  <  Oq  <  6^  -  <g,  we  obtain 


lim  ^oY('‘o) 

V  =  .Y-0'^ - 71 


-Sg  (^2  +  52  ^  ^2  ,2^  -  I 


lim 

X-*Q 


V  - 


4170"^ 


dS 


-j 


-Sq  [XU  5^  +  a"  (0/-  V  I 


2,1/2 


d8 


-Sq  [X^^S^^aU0i-e^-2»^f] 


1/2 


1 1  (»b«  ^0>  =  ®  ^0  =  -  ‘o  * 

Integrating  and  taking  the  limit,  we  finally  obtain 


or 


'"[  *0 * ^ (»| - "o)'  ] 

-In  [-S,  +  +  n^  (9,  -  j 

-ln[  S„  +  y«o' +  o“  («,  -  »o  -  2  ] 

+  ln  [rS(,  +  \/«o  +  (»,  -  -  2  *1))^  j| 

'o>'<'o)  r  /-^ - 5 - jI 

I"[-«<,  +  y«o  +"  («|-V 


^1  (’o’  ^0^ 


since  the  second  term  is  certainly  the  most  important  term  as  9^  -*  d^.  As  expected,  the 
induced  velocity  components  have  logarithmic  singularity  at  the  leading  edge. 

Equation  [16b] 

^bt  (’o’  (’b’ 


lim 

=  X-OV 


-’b  +  ^o  y(r)i9-9^)lpcosid-9^)  +  p'(r^-rco3{e-9^))] 


^11 


‘^'b“^o  ^o~‘o  +  cos(0-0o)  +  (p/;7)2(0-0p)2] 

Following  the  same  procedure  and  neglecting  the  p'torm,  we  obtain 

Wo)  -  -  1"  F-  *0  *  y^o 

40"'  rz  ^  J 


3/2 


do  dr 
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Equation  [16c] 


^/a(V^o)  =  '3(VV 

lim  f  1  A""*®  r^(r,9)f  Ir-fjj  cos  (0-0q)1 

J  J - 


1  'b~^o  ^o~*o  W^  +  »-o^  +  '-^-2vcos(0-0q)  +  (p/;7)2(0-0q)2]^'^ 


reir,d)  =  rQ(r^,d)^r0ir^,e){r-r^) 


we  have 


lim 

^3  ^*'o’  ^0^  ~  i^-*0 


r^o  /.^o 


^  r  r  °  rg(r^,o..rg-(r,,.)g 

4n  J  J  T  o  -/  3/2 

"o  -S„  MX’ +  8= + 


5  dddt 


r  'o  f‘°  f*“  r'^Vo.‘)  +  C(r„)8-y'(r„).S 


lin,  I  '0  f  °  f 

x-oi-JT  J  J 

I  *0  "^0 


iX^  .  8^  . 


8  d8di 


since  r^‘'(^'o>*)  can  be  expressed  as  the  sum  of  C{rQ)  and  -y'{rQ)*  where  C’(rQ)  represents 
quantities  evaluated  at  r^,  a  constant. 

The  integration  of  the  first  term  with  respect  to  8  and  that  of  the  third  term  with  respect 
to  i  vanish,  and  the  only  remaining  term  is  the  second  term.  Therefore 


lim  f  f‘” 

'sCo'V'X-O-l  —  J  J. 


d8dt 


-4^  -5,  iX^^8UaU^) 
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Her-^  again  we  observe  that  the  most  important  term  is  the  last  term.  Carrying  out  the 
integration,  taking  the  limit,  and  using  the  approximation  of 


-S, 


+  (S  2  + 


1/2 


we  finally  have 


^3  (^0>  ^o)  = 


4rr 


“1 


<  In  -  +  24  (In  4e  -  1) 

2So 


which  vanishes  as  Cq  -♦  0. 


Equation  [16cl] 


lim 


’b'^'^o  [f^_rcos(0-0Q)-f(0-0o)sin(0-0p)] 


^  -  V  up  V  pr^(rQ,0)ir^-rc 

U^+,„ni-^-2r„r 


3/2 


r 


lim 

'■x-*o 


*0  r^O 


r^rco3{e-e^)+(]x/ny{e-0Qy]  J 


yd$dr 


-f  I 


“•o  "“^0 


d8di 


which  vanishes  by  following  the  same  produced  as  for  ly 

This  shows  that  the  only  contribution  is  from  /j  and  /j. 


CONTRIBUTION  TO  MEAN  LINE 

The  contribution  to  the  mean  line  due  to  the  small  area  in  the  close  neighborhood  of 
the  singular  point  when  0i>  0q>  6^  -  is 


D  r„ 


*0) 


2V  cos  /3,. 


f'' 
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where  is  the  induced  downwash,  which  can  be  expressed  as 


('■o’  ^0^  =  ^  (‘‘o’  ^o)  ^.('■o)  -  ^2  ('‘o’  ^o) 


Substituting  the  results  of  /j  and  and  making  the  following  approximation 


-5oW5o'  +  «'  -  So  — 


2  \  5, 


we  get 


^  ('■o’  ^o)  =  ^('■o) 


2  \  5. 


where 


^(*■0)  2 

/('■q)  =  — [»'o  ^i('’o)  +  (P('‘o>/'^)  /3,-(»‘o)l 

4  a  77 


Performing  the  integration,  we  obtain 


Z?  f. 


^^'"('■o’V=  2VCOS/3. 


■  ^o«'  /  *0 

In - +  2  (  In - i 

2  \  5o  - 


Since  <<  1,  this- contribution  to  the  mean  line  is  indeed  negligible.  An  estimate 
for  a  case  where  =  0.001  and  5^  =  0.005  shows  that  the  contribution  of  this  small  area  as 
compared  to  that  of  the  remaining  area  is  in  the  order  of  tenths  of  1  percent.  These  values  of 
^0  and  §0  actually  used  in  the  computing  program. 
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APPENDIX  B 


PREPARATION  OF  INPUT  DATA 


The  input  data  consists  mainly  of  two  parts:  (1)  program  controls  and  (2)  propeller 
data.  The  following  discussion  is  outlined  in  the  order  of  the  input  data  format;  see  Figure 
16. 


Card  1  Format  918 

NL  Number  of  radial  stations  in  the  common  network  of  field  points  (maximum  number 

of  radial  stations  in  the  tip  or  hub  zones) 

Specify  NL  =  18 

Card  2  Format  8F9.5 

DY  Sq,  half-width  of  mid-zone 

Specify  DY  =  0.005 

DFD  Minimum  chordwise  spacing,  in  the  common  network  of  field  points  in 

degrees 

Specify  DFD  =  0.1* 

FDFD  Maximum  chordwise  spacing,  in  the  common  network  of  field  points  in 

degrees 

Specify  FDFD  =  2.0* 

DFFD  Maximum  range  of  <f>  in  the  common  network  of  field  points  in  degrees 
Specify  DFFD  =  150.0* 

Card  3  (more  than  one  card)  Format  9F8.5 

SPR  Spacings  of  radial  stations  in  the  common  network  of  field  points.  Total  number 
of  entries  should  equal  NL,  The  following  set  of  numbers  may  be  specified: 
0.005,  0.010,  0,016,  0.024,  0.C35,  0.050,  0.070,  0.096,  0.129,  0.170,  0.220,  0.280, 
0.352,  0.440,  0.552,  0.709,  0.945,  1.000 

If  Cards  1,  2,  and  3  as  specified  are  used,  no  entry  for  these  cards  is  necessary. 

Card  4  Format  12A6 

Propeller  identification  not  more  than  72  digits. 

(This  is  the  first  line  in  the  print  out.) 

Card  5  Format  9l8 

NN  Number  of  blades 

NKN  Number  of  chordwise  coordinates  at  which  induced  velocities  are  calculated. 


*Total  number  of  chordwise  stations  on  either  side  of  control  point  should  not  exceed  90, 
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NCL  Number  of  chordwise  coordinates  x/Lj.  at  which  camber  distribution  is  calculated. 
NCTR  Program  control. 

When  calculation  for  a  given  set  of  input  data  is  completed,  if: 

NCTR  =  1,  the  computer  reads  another  completed  set  of  input  data  (cards  4 
through  9)  and  proceeds  to  compute. 

NCTR  =-  1,  the  computer  reads  another  set  of  data  for  NN,  NCTR,  NPT,  KSP, 

VM,  DY,  SM,  SW  and  proceeds  to  compute. 

NCTR  =  0,  program  terminates. 

KRN  Number  of  radial  coordinates  at  which  induced  velocities  and  camber  distribution 
are  calculated. 

NPT  Specify  NPT  =  0.  Program  control  for  print  out.  (See  Appendix  C.) 

KSP  Specify  KSP  =  0.  Program  conti'ol  for  print  out.  (See  Appendix  C.) 

Card  6  Format  9F8.5 

DK  Used  to  calculate  chordwise  coordinates  at  which  induced  velocities  are  calcu¬ 

lated.  The  chordwise  coordinate  is  calculated  by  using  the  formula: 

d/  =  d^(g-DK(;) 

Total  number  of  entries  should  equal  NKN. 


Card  7  Format  9F8.5 

XL  Chordwise  coordinates  at  which  camber  distribution  is  calculated. 

Total  number  of  entries  should  equal  NCL. 

Card  8  (9  cards)  Format  9F8.5 

Propeller  data: 

1  st  card  -  Radii  r  at  which  input  data  is  given  in  the  following  four  cards. 

2nd  card  -  Dimensionless  leading  edge  coordinate  from  reference  line  at  r,  L^/D. 

3rd  card  -  Dimensionless  chord  length  at  r,  Lj,/D. 

4  th  card  —  Tangent  of  hydrodynamic  pitch  angle  at  r,  tan  taken  from  ohe  lifting  line 

calculation  output.* 

5  th  card  -  Dimensionless  bound  circulation  distribution  r(r)  based  on  ship  speed  taken 

from  the  lifting-line  calculation  output.* 

6  th  card  -  Radii  at  which  induced  velocities  and  camber  distribution  are  calculated. 

Total  number  of  entries  should  equal  KRN. 

7th  card  -  Advance  velocity  at  r^,  V 

8  th  card  -  Axial  induced  velocity  at  taken  from  the  lifting-line  calculation  output.* 

9  th  card  -  Tangential  induced  velocity  at  taken  from  the  lifting-line  calculation  output.* 


*AML  Problem  No.  106,  "Moderately  Loaded  FVopellers  Using  Induction  Factors." 
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Card  9  Format  9F8.5 


VM  Scale  factor  for  chord  length  used  for  calculating  the  camber  for  different  expanded 

area  ratios  with  respect  to  a  reference  EAR. 

Specify  VM  =  1.0  if  the  input  data  Lj./D  (third  card  of  Card  8)  give  the  required 
EAR.  For  other  values  of  EAR  it  is  only  necessary  to  change  the  value  of  VM, 
e.g.,  VM  =  1.25  represents  an  EAR  which  is  25  percent  larger. 

AJ  Coefficient  of  advance,  based  on  ship  velocity. 

DiA  Propeller  diameter. 

SM  Specify  SM  =  1.0. 

If  only  the  camber  due  to  bound  circulation  is  to  appear  in  the  print  out, 
specify  SM  =  0. 

SW  Specify  SW  =  1.0. 

If  only  the  camber  due  to  free  circulation  is  to  appear  in  the  print  out, 
specify  SW  =  0. 

Note  1.  The  program  is  written  so  that  more  than  one  problem  can  be  handled  at  one  time. 

The  computer  performs  complete  sets  of  calculations,  one  for  each  problem, 
successively.  The  problem  could  be 

a.  for  entirely  different  propellers,  in  which  case  one  set  of  the  input 
data  Cards  4  through  9  for  each  propeller  v  ill  be  required, 

b.  for  a  propeller  the  same  as  the  one  specified  in  the  immediate 
preceding  set  except  that  the  expanded  area  ratio  (EAR)  is 
different, 

c.  for  the  same  propeller  as  specified  in  the  immediate  preceding  set 
except  calculation  is  to  be  performed  with  different  NN,  NPT,  KSP 
or  SM,  or 

d.  for  a  combination  of  b  and  c. 

For  either  b,  c,  or  d,  a  set  of  two  input  cards  will  be  required.  They  specify  the 
quantities  involved  in  the  following  order: 

Card  i  NN,  NCTR,  NPT,  KSP  (Format  418) 

Card  ii  VM,  DY,  SM,  SW  (Format  4F8.5). 

In  each  case,  apfx’opriate  values  of  NCTR  should  be  specified.  The  sets  of  input 
cards  are  to  be  arranged  consecutively  in  an  appropriate  order. 

Note  2.  The  computing  time  depends  mainly  on  the  number  of  blades,  the  shape  of  blade 
outline,  and  the  fineness  of  network  of  field  points,  i.e.,  the  number  of  stations. 

For  a  conventional  propeller  as  the  one  shown  in  Figure  16,  the  computing  time 
required  is  tabulated  for  various  values  for  Z,  NL,  and  maximum 

range  of  4>' 


z 

NL. 

(A^) 

'  ^'max 

Range  of  (f> 

Time 

1 

20 

0.10° 

2.0° 

150° 

9.4  min 

1 

11 

0.10° 

2.0° 

150° 

6.0  min 

1 

11 

0.10° 

3.0° 

220° 

5.5  min 

5 

20 

0.10° 

2.0° 

150° 

11.4  min 

46 


9«t«* 

flMT  OATA 


t« 

•••0900  O.IOOOO  2*00000190. 00090 


••00900 

0.01000 

o.oieoo 

0.02400 

0.03900 

0.09000 

0.07000 

••17000 

0*22009 

0.20000 

0.39200 

0.44000 

0.99200 

0.70900 

9 

• 

o 

0 

o 

6 

3 

••OOffS 

0.090«S 

0.11062 

0.20414 

0.20966 

0.301 ja 

0.449t7 

••oe«ts 

0.01290 

0.02990 

0.09000 

0.10000 

0.20000 

0.30000 

••20000 

0.30000 

0.40990 

0.90000 

0.60000 

4.70000 

o.aoooo 

••ISAOO 

0. 10600 

0.21990 

0.23100 

0*23200 

0.21790 

0.10900 

••30000 

0.37200 

0.43000 

0.46300 

0.46400 

0.43900 

0.37000 

i*»42ae 

I. 32039 

0.00766 

0.70420 

0.96622 

0.46392 

0.30696 

•  • 

••01730 

0.02266 

0.02420 

0.02203 

0.01936 

0.01463 

••29000 

0*30900 

0.40000 

0.90000 

O.6CC00 

0.70000 

o.eoeoo 

l•60000 

t.ooooo 

1.09000 

1.00000 

1.00000 

1.00000 

1.00000 

••14900 

O.I6902 

0.20661 

0.22173 

0.20769 

0.17661 

0.13904 

••16010 

e«t6670 

0.16039 

0.13704 

0.10917 

0.07403 

0.04089 

t^MAOO 

•••2600 

0.00000 

i.eoooo 

1.00000 

(•*c  Appendix  9  for  czpt«estlo'> 


o.ooaoo-o« 

O.OCOOO  0.09000 

t.ooooo  t.ooooo 
o.oaosr  v.oaooo 

0.02767  0.01770 

Card  9 


AO^ieO  MATHCMATKS  tAaOAATOA  XOLU 

mO^CLLfO  0C9ICN  OASCD  On  Cl^tlNC  SUA^ACC  T*^0AT 
VNtrOftH  ChOAObISC  LOaO  OIStCieUT tON 

AAI.U  Card  4 


Card 

1 

Card 

2 

0.09600 

0.12900 

Card 

3 

0.94900 

1.00000 

0*49007 

Card 

} 

0.40000 

0.90000 

Card 

6 

0.90000 

1.00000 

Card 

7 

0.13400 

0.02000 

Card 

S 

0.26000 

0.04000 

O. 32722 

0.27960 

r 

L/Cr)/I> 
lt(r)/D 
tan  ^i(r) 

r(r) 

v,(r.) 

utfro) 


•wmiT 

•  XIAL  tmvCMO  VCLOCITT  CO.MM’HT 


•  X/AO 

0.2900 

0.3600 

0.4000 

0.9600 

0.6000 

0.7000 

0.0000 

0.9000 

0.9900 

*✓1.7 

e.0999 

0.0990 

0«I047 

e.i3i4 

0.1941 

0.1922 

0.1460 

0.  1389 

0.1299 

0.1146 

0.0900 

0.0421 

0.0769 

0.0903 

0*1091 

0.1064 

0*0999 

0.0994 

0.0007 

0.0919 

0*tl06 

0.0309 

0.0960 

0.0662 

0.07&9 

0.0774 

0.0662 

0.0691 

0.0665 

0.0724 

0.2041 

0.0214 

0.0403 

0.0440 

0*0940 

0.0924 

0.0400 

0.0466 

0.0461 

0.0970 

0.2999 

0.0141 

0.0264 

0.0209 

0.0343 

0.0337 

0.0313 

0.0303 

0.0324 

0.0372 

0.3014 

0.0003 

0.0149 

0.0194 

0.0191 

0.0166 

0.0179 

0.0172 

0.0166 

0.C249 

••4492 

0.0043 

0.0064 

0.0061 

0*0902 

0.0079 

0.0074 

0*0074 

0.0079 

0.0109 

0.4901 

0.0019 

0.0013 

0.0006 

0.0010 

0.0016 

0.0019 

0.0019 

0.0013 

0.0026 

TAMOCMTtAt.  iNOOCeC  VCLCCITY 

C0MAOt<NT 

AX/AO 

0.2906 

0.3000 

0*4000 

0.9COO 

0.6000 

0.7000 

0.0000 

0.9000 

0.9900 

A/UT 

0«0099 

-0.1699 

-0.1993 

-0*1730 

-0.1400 

-0*1142 

-0*0049 

-0.0620 

-0.049.’ 

•0.0360 

0.0900 

-0.1404 

-0.1970 

-0*1393 

-0.1147 

•0.0677 

-0.0623 

-0.0447 

•0.0333 

-0.0290 

0*1 106 

-0.1190 

-0.1277 

-0*1073 

-0.0009 

-0.0671 

-0.0494 

•0.0323 

-0.0293 

-0.0230 

0.2041 

-0.0009 

-0.0974 

-0*0809 

-0.0664 

-0.C466 

-0.0336 

-0.0232 

-0.0104 

-0.0100 

0«2999 

-0.0614 

-0.0671 

-0.0947 

-0.0447 

-0.C323 

-0.0220 

-0.0193 

-0.0124 

-0.01 19 

0.3614 

-0.0363 

-0.0,391 

-0.0319 

-0.0257 

-0.0104 

-0.0129 

•0.0007 

-0.0071 

-0.0070 

0.4492 

-0.0163 

-0.0160 

-0.0132 

-O.OItt 

-0.0070 

-0.0993 

-0.0030 

-10.0030 

-0.0039 

0.4901 

-0.0043 

-0.0034 

-0.0022 

-0.0024 

-0.0815 

-0.0010 

-0.0000 

-0.0009 

-0.0009 

fMOUCCO  OOVMVASH  CVC  TO  Al.i. 

8CA0e9 

AX/AO 

0.2900 

0.3000 

0.4000 

0.9COO 

0.9000 

0.7000 

0.0000 

0.9000 

0.9900 

X/CT 

••••99 

0.1740 

0.2.179 

0.2136 

0.2112 

0.1007 

0.1602 

0.I9IA 

0.1339 

0.1201 

••% v09 

0.1400 

0.'l707 

0.1970 

9.1993 

0.1376 

0.1164 

0.1091 

0.0947 

0.0964 

••tlAO 

0*1 129 

0.1340 

0.1212 

0.1132 

0.1004 

0.0791 

0.9724 

0.0710 

0.0799 

••2041 

0.0094 

0.1010 

C.0073 

0.0824 

0.0696 

0.0904 

0.0919 

O.09I4 

0.0990 

9.2999 

c.o9oe 

0.0607 

0.0979 

0.093a 

0.0493 

0.0376 

0.0337 

0.0347 

C.039I 

0. 3014 

0.0340 

9.0397 

0.0329 

0.0304 

0.0294 

0.0211 

0.0192 

0.0199 

0.0297 

•  .44  72 

0.0159 

0.01  \ 

0.0134 

0.0131 

0.0107 

0.0090 

0.0062 

0.0064 

0.0119 

••4901 

0.0046 

O.CV.«9 

0.0020 

0*0029 

0.0021 

0.0016 

0.001 7 

0.0014 

0.0027 

AX/AO 

0.2900 

0.3000 

0.4000 

0.9''00 

0.6000 

0.7000 

0.0000 

0*9000 

0.9900 

A/0 

1.1960 

1.1977 

1.1406 

1.1063 

1.0073 

1.0202 

0.9726 

0.9292 

0.0976 

L.d. 

1.3933 

1.4600 

1.7240 

1.6480 

1.0960 

1*7400 

1  .4000 

1 .0720 

0.7492 

CHOAO 

2.7213 

2.9760 

3.4400 

3.9960 

3.7120 

3.4000 

2.9600 

2.1440 

1.4943 

CAAtCA  OtSTAieuTION 

X/CT 

••0099 

0.0039 

0.0044 

0.0042 

0.0036 

0.0029 

0.0022 

0.0019 

0.0006 

0.0009 

0.0129 

0.0044 

0.0095 

0.0052 

0.0047 

0.0037 

0*0027 

0.0016 

O.OOt 1 

0.0006 

0.9290 

0.0004 

0.0109 

0.0099 

0.0069 

0.0069 

0.0090 

0.0034 

0.0020 

0.0012 

0.0900 

0.0197 

0.0194 

o.oiao 

0.0162 

0.0129 

0.0090 

0.0061 

0.0039 

0.0022 

0.1000 

0.0203 

0.0346 

0.0314 

0.0260 

0.0216 

0.0191 

0.0103 

0.0061 

0.0039 

0.2000 

0.0403 

0.0909 

0.0919 

0.0499 

0.0340 

0.0238 

0.0163 

0.0099 

0.0067 

0.3000 

0.0624 

0.0791 

0.0696 

0.0572 

9.0A33 

0.0297 

0.0202 

0.0129 

0.0007 

0.4000 

0.0709 

0.0049 

0.0739 

0.0639 

0.0402 

0.0330 

0.0229 

0.0140 

0.0099 

••9090 

0.0740 

0.0002 

0.0760 

0.0661 

0.0497 

0.0341 

0.0232 

0.0149 

0*0104 

CAA«CA  AATIO 

X/LT 

0.0009 

0.0013 

0.0019 

0.0012 

0.0010 

0.0000 

0.0006 

0.0005 

0.0004 

0.0003 

0.0129 

0.0016 

0.0010 

0.0019 

0.0013 

0.0010 

0.0000 

0.0006 

o.ooos 

0.0004 

0.0290 

0.0031 

0.0039 

0.0029 

0.0024 

0.0019 

0.0014 

0.0012 

0.0009 

0.0000 

0.0900 

0.0090 

0.0069 

0.0052 

0.0944 

0.0034 

0.0026 

0.0021 

0.0016 

0.0019 

O.IOOO 

0.0104 

0.0116 

0.0091 

0.0076 

0.0090 

0.0043 

0.0039 

0.002C 

0.0026 

0.2000 

0*0177 

0.0197 

C.019I 

0.0123 

0.0094 

0.0060 

0.0099 

0.0046 

0*0049 

O.SOOC 

0.0229 

0.0292 

9.0191 

0.0195 

0.0117 

0.0009 

0 .0060 

0.0098 

0.0090 

9.4000 

0.0261 

0.0209 

0.0214 

0.0173 

0.0130 

0.0099 

0.0076 

0.0069 

0.0066 

•••200 

0.0272 

0.0296 

0.0221 

0.0179 

0.0134 

0.0090 

0.0070 

0.0066 

0.0069 

Figure  16  -  Sample  Computer  Output  (for  Propeller  3916A) 
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APPENDIX  C 


SAMPLES  OF  OUTPUT 


The  program  prints  out  the  following  information  in  tabular  form;  see  Figure  16. 

Axial  induced  velocity  component. 

Tangential  induced  velocity  component. 

Induced  downwash. 

Camber  distribution  (in  unit  same  as  propeller  diameter). 

Camber  ratio. 


In  addition,  the  program  has  a  feature  of  printing  out  certain  intermediate  information 
when  desired;  see  Figure  17.  It  will  print  out  the  chordwise  integrated  values  of  the  axial 
and  tangential  induced  velocity  components  in  the  following  general  order: 


Tip-zone 

r  =  ri 

chordwise  integrated  values  of  at  various  6^ 
chordwise  integrated  values  of  at  various 


(same  print  out  as  for  f  =  fj) 


etc. 

Hub-zone 

(same  print  out  as  in  the  tip-zone) 
etc. 

Integrated  values  of  in  mid-zone  at  various  Oq 
Integrated  values  of  in  mid-zone  at  various  6^ 

At  r=f2 

(same  print  out  for  Tq  =  ) 


etc. 


Also  printed  out  are  the  values  of  integrands  in  the  equations  expressing  the  induced 
velocity  components  (Equation  [16])  for  any  specified  control  point  Pir*,  0^)  in  the  following 
order: 
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At  r  =  and  6^  =  0^ 

Tip-  and  Hub-Zones 
At  r  =  Tj 

Chordwise  stations  6,  PHI 

Axial  components  due  to  superimposed  lifting  line  (to  be  exact,  the  second  term 
of  Equation  [16c])  corresponding  to  UHH 

Tangential  components  due  to  superimposed  lifting  line  (the  st-cond  term  of  Equa¬ 
tion  [16d])  corresponding  to  6,  VHH 

Axial  components  due  to  free  system  (the  first  term  of  Equation  [16c])  corresponding 
to  ^  —  leading  edge  portion  (from  P  to  leading  edge),  UHL 

Axial  component  due  to  free  system  —  trailing  edge  portion  (from  P  to  trailing 
edge),  UHT 

Tangential  component  due  to  free  system  (the  first  term  of  Equation  [16d]) 
corresponding  to  ^  —  leading  edge  portion,  VHL 

Tangeniial  component  due  to  free  system  —  trailing  edge  portion,  VHT 

Axial  component  due  to  bound  system  (Equation  [16a])  corresponding  to  — 
leading  edge  portion,  UBL 

Axial  component  due  to  bound  system  —  trailing  edge  portion,  UBT 

Tangential  component  due  to  bound  system  (Equation  [16b])  corresponding  to 
<f>  8c  leading  edge  portion,  VBL 

Tangential  component  due  to  bound  system  —  trailing  edge  portion,  VBT 

Value  of  f  j 

Chordwise  integrated  values  of  at  various  Oq 

Chordwise  integrated  values  of  v^&t  various  6^ 


At  r  =  r2 

(similar  print  out  as  for  r  =  fj) 


Mid-Zone 

(similar  print  out  for  the  axial  and  tangential  components  as  in  the  tip-  and  hub-zones. 
Integrated  values  of  various  Oq 
Integrated  values  of  at  various  6q 

Figure  16  shows  the  output  for  a  typical  problem. 

Figure  17  shows  the  partial  output  of  the  same  problem  as  in  Figure  16,  except  inter¬ 
mediate  information  for  point  P  at  =  r®  =  0.7  and  6^  =  6^,  i.e.,  NPT  =  6,  KSP  =  3,  is  also 
printed. 
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CHOHOalSe  INTE6KATE0  VALUES  OF  VA  AHO  VT  AT  RO'>.2SaO 


TIF-ZONE 

Ra.tSSO 


02 

-0*889206 

C2 

-0*762316 

02 

-0*592216 

02 

0*l24*3e 

03 

0*120396 

03 

0*103106 

03 

0*799496 

02 

R**2ftU0 

-0.43474e 

02 

-0*424196 

02 

-0*367326 

02 

-0*286356 

02 

e«s349ie 

02 

0*926746 

02 

0*495976 

02 

0.35442E 

02 

*0*274S0e 

02 

-0*269746 

02 

-0*239706 

02 

-0*184406 

02 

02 

0*30613E 

02 

0.26S036 

02 

0*2089 76 

02 

-0«t«773e 

02 

-0*182446 

02 

-0*160816 

02 

-0.126346 

02 

0«>««13C 

02 

0*184776 

02 

0*164046 

02 

0*128486 

02 

INi«2790 

-0*12769C 

02 

-0*123486 

02 

-0*109766 

02 

-0*866976 

01 

0.10670E 

02 

0*1062ftE 

02 

0*960996 

01 

0*758646 

01 

K««2900 

-0*9«030E 

Ot 

-0*829476 

01 

-0*738986 

01 

-0*587246 

01 

••SSS4CC 

01 

0*994386 

01 

0*514226 

01 

0*41*1756 

01 

R«*3C30 

-04.57757e 

01 

-0*590016 

01 

-0*493646 

01 

-0.395436 

Ot 

e»746S4E 

01 

0*248046 

01 

0*239456 

01 

0*197106 

oi 

n•^'S^50 

-0.3\34SC 

01 

-0*372606 

01 

-0*334436 

01 

-0*270026 

01 

Q*7d56SC 

00 

0*818506 

00 

0*866676 

00 

0*762496 

00 

IN)«3S10 

-0.27SS3C 

01 

-0.360696 

01 

-0*233986 

01 

-0*190636 

01 

-0*2}427E- 

l»«*3«40 

-01 

0*195076 

-01 

0*109466* 

-00 

0*1494  76- 

-00 

->0*t94«6E 

01 

-0*184396 

01 

-0*166126 

01 

-0*136996 

01 

-o*3«ta9e-oo 

-0.32030E- 

-09 

-0*227706-00 

-0*135626-00 

fl-•4^30 

-0*t2396C 

01 

-0*118256 

Oa 

-0*107766 

01 

-0*905686 

00 

-0«4«S37E- 

ll■«47S0 

-00 

-0*448286- 

-QU 

-0*366076* 

-00 

-0*262836- 

-00 

-9«4«479E' 

-00 

-0*474376 

-00 

-0*451436* 

-00 

-0*399826-00 

-0*309342 

00 

-0*478526 

-00 

-0*410426-00 

-0*313076-00 

Ma^SSSO 

0*24920E-00 

0*246226- 

-00 

0*187726-00 

0*121586-00 

-0*47020e‘ 

-00 

-0*448326-00 

-0.394216* 

-00 

-0*310596- 

-00 

«a*6070 

0*077I5£ 

00 

0*803786 

00 

0*692606 

00 

0*532076 

00 

-0*300S0E' 

-00 

-0*368286 

-00 

-0.328796-00 

-0*262986-00 

9«*4990 

e*t2sai€ 

01 

0*119886 

01 

0*977086 

00 

0*782586 

00 

-0*278722-00 

-0*263276-00 

-0*231166-00 

-0*187906-00 

8««8070 

0«S4985E 

01 

0*137916 

01 

0*117876 

01 

0*919456 

00 

-0«tA82Se-00 

-0*197806-00 

-0*138446-00 

-0*110826-00 

8««9840 

0*22769E 

01 

0*212936 

01 

0*181866 

01 

0*141916 

01 

-0*96S2tE-0t 

-0*901026-01 

-0.787286-01 

-0*631176-01 

na*9999 

0*39923e 

01 

0*326486 

01 

0*278796 

01 

0*217306 

01 

-0*12182E-00 

-0*113586-00 

-0*992796-01 

-0*795266-01 

HilO-ZOHE 

8««2490 

0«tll91E 

03 

0*967706 

02 

0*812726 

02 

0*627406 

02 

-0*177a2C 

aa*2400 

03 

-0*154636 

03 

-0* 130216 

03 

-0.100756 

03 

o*saa9iE 

02 

0*504876 

02 

0*420746 

02 

0*323966 

02 

-o.iooaiE 

Ra«2390 

03 

-0*874876 

02 

-0.733736 

02 

-0.567696 

02 

0«39929E 

02 

0*350006 

02 

0*290406 

02 

0*223176 

02 

-0*74085E 

a«*2290 

02 

-0*657686 

02 

-0*551286 

02 

-0*426846 

02 

0*2a9a8E 

02 

0*260336 

02 

0*215976 

02 

0.165746 

02 

-0*S9481E 

fH>*22tO 

02 

-0*539986 

02 

-0*454226 

02 

-0*352316 

02 

0*21292E 

02 

0*196716 

02 

0*164406 

02 

0*126186 

02 

-0«S0409E 

R«*2I00 

02 

-0*467916 

02 

-0*397266 

02 

-0*309216 

02 

0*19763E 

02 

0*  148946 

02 

0*126406 

02 

0.973O46 

01 

-0*499776 

aa*2001 

02 

-0*432126 

02 

-0*372206 

02 

-0*291406 

02 

0«t2794E 

02 

0*121896 

02 

0*105586 

02 

0*810206 

Ot 

-0*49627E 

02 

-0*430756 

02 

-0*374906 

02 

-0.294956 

02 

IHTECRATEO  VA.VT  IN  MIO-ZONE 
0.43029E-02  -0.780106-03  -0.128286-02  -0.114506-02 
-0.82402E-02  -0.593326-03  0.411776-03  0.343126-03 


-0« 407146  02 
0*548846  02 

-0*233736  02 
0*314796  02 

-0*960426  01 
0.129266  02 

-0*129346  Ot 

0*173366  01 

-0*196906  02 
0*243036  02 

-0.112646  02 
0*135776  02 

-0*456216  01 
0*561116  01 

-0*515306  00 
0*623836  00 

-0*126816  02 
0*143176  02 

-0.722846  01 
0*813786  01 

-0*208396  01 
0*323616  01 

-0.2S932E-00 

0*278706-00 

-0*869136  01 
0*879646  01 

-0*493276  01 
0*497226  01 

-0*193216  01 
0*193566  01 

-0*118516-00 
0* 103246-00 

-0*596736  01 
0*519296  01 

-0*3^.6876  Ol 
0*291386  01 

-0*128856  01 

0*110086  01 

-0*304046-01 

0*643866-02 

-0*404706  01 
0*282266  01 

-0.227156  01 
0*156896  01 

-0.344866  00 
0*568076  00 

0*166076-01 

-0*368906-01 

-0*273336  01 
0.135996  01 

-0*152896  01 
0*748066  00 

-0*556466  00 
0*256146-00 

0.327996-01 

-0*414026-01 

-0*187946  01 
0*537966  00 

-0.105676  01 
0*294086-00 

-0*388386-00 

0*949126-01 

0.171916-01 

-0*261696-01 

-0*134456  01 
0*1200?6-00 

-0*771536  00 
0*670506-01 

-0*302336-00 

0*209696-01 

-0*168206-01 

-0*709126-02 

-0*985606  00 
-0*775936-01 

-0*585116  00 
-0.403236-01 

-0*253166-00 

-0*131056-01 

-0*503106-01 

0*302706-02 

-0*668416  00 
-0*170136-00 

-0*411636-00 
-0. 932146-01 

-0*194346-00 

-0.333226-01 

-0*605736-01 

0*297256-02 

-0*309656-00 

-0*210856-00 

-0*198716-00 
-0* 1 19156-00 

-0*992836-01 

-0*459866-01 

-0*369326-01 

-0.335926-02 

0*672306-01 

-0*213816-00 

0*320206-01 

-0*122916-00 

0*114636-01 

-0*512126-01 

0*620076-03 

-0*804206-02 

0*361936-00 

-0*182926-00 

0*207956-00 

-0*105516-00 

0*883756-01 

-0*440636-01 

0*175086-01 

-0*738626-02 

0*539126  00 
-0*131596-00 

0*309836-00 

-0*761426-01 

0.127716-00 

-0*316596-01 

0*178306-01 

-0*480096-02 

0*649756  00 
-0.793736-01 

0*377496-00 

-0*466706-01 

0* A 59606-00 
-0*197406-01 

0*282096-01 

-0*334926-02 

0*979156  00 
-0*446696-01 

0*561036  00 
-0*260206-01 

0*242086-00 
-0*1 12846-01 

0*408196-01 
-0. 190046-02 

0*150506  OS 
-0*564076-01 

0*887926  00 
-0.337896-01 

0*382836-00 
-0* 146646-01 

0*575016-01 

-0*220536-02 

0*431196  02 
-0*693546  02 

0*249106  02 
-0*401036  02 

0*104996  02 
-0*169106  02 

0*100786  01 
-0*291496  01 

0*222616  02 
-0*391416  02 

0*120936  02 
-0*227126  02 

0*549026  01 
-0*967956  01 

0*102786  01 
-0*181366  Ot 

0.153316  02 
-0.294786  02 

0*089746  01 
-0.171576  02 

0*362046  01 
-0.737506  01 

0*761146  00 
-0*146976  01 

0*113836  02 
-0.-243776  02 

0*661666  01 
-0*142206  02 

0*286076  01 
-0.616076  01 

0*598526  00 
-0*128056  01 

0.C6615E  01 
-0*214436  02 

0*503656  01 
-0*125406  02 

0*216416  01 
-0*544996  01 

0*467016-00 
-0*116456  01 

0*667626  01 
-0.202476  02 

0*386006  01 
-0*110236  02 

0*166026  01 
-0*509036  01 

0*331666-00 
-0*101696  01 

0*560046  01 
-0.204876  02 

0*321606  01 
-0*118786  02 

0*134166  0! 
-0.49729E  01 

0*213296-00 
-0*791246  00 

-0*829646-03 

0*436646-03 

-0*467876-03 

0*266576-03 

-0*204016-03 

0*112506-03 

-0*305326-04 

0*162356-04 

(Similar  print  out  for  Ro»0.3,  0.4,  0.9S) 

Figure  17  -  Sample  Computer  Output  with  Intermediate  Information 
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CMOHOSfM  SNTCQIATCO  VAtUCS  VA  AND  VT  AT  ftO«*7QOO 

YI^ZOMC 

THCTASUao  3 

iMTeimeOIATe  XHTOAMATIOM  of  PMl«UHH»VKH.U«.*OHT*VHL*VHT«Oet.«U«T»va(.*V8T  AT  A  t 
•«  0«1745AC-02  0*5236ll'-'02  0*I0472C-0|  0*X74SAC-0t  0«2618dC«0t  0*38652e-0l  0«4eB70e'0t  O.62ei3C'0t 

••76S4ie*0l  0«^e«4e*0l  O.tlSl^e-OO  0*t2762C-00  e«l^Be3e-00  0«te32frC-00  O«20944£>00  0.2j737e>00  0*26704C>00 

t*2«SAAe*00  0«33t62C-00  0.3AB92C-00  0«40I43€-00  0«43634C>00  0*47t2S£-00  O.SOfrtSE  00  0.54I0BC  00  0.57507C  00 

••OtOOTC  00  0«6AS7«C  00  0.6«069e  00  0.7I35OC  00  0*7SOSOE  00  0*7a5AtC  00  0«52032r  00  0«aS522e  00  o.aooue  00 

••OZSOAC  00  0.9S094e  00  0.90762e  00 

•'•749072  04  0.674932  04  0«3S233e  04  O.IOTOIC  04  0*32t20C  03  ••10672C  03  0.43A90C  02  0.20254C  02  O.lOTazC  02 

••643132  01  0.426022  Ot  0^30716C  OI  0.260S7C  01  0.192442  Ot  0.163132  01  0.t42a3£  01  0.128162  01  0.I1723C  01 

••10884C  Ot  0.102332  01  0.970942  00  0.930232  00  0.697742  00  0.a7109E  00  0.646672  00  0.829262  00  0.811682  00 

••793712  00  0.7t00S2  OO'  0.764242  00  0.747692  00  0.729892  00  0.710392  00  0.686922  00  0.663362  00  0.639632  00 

9.612772  00  0.384802  00  0.362332  00 

-0.344442  04  -0.310092  04  -0.161862  04  -0.302732  03  -0.147142  03  -0.494922  02  -0.193132  02  -0.683412  01  -0.447862  01 

-e.248S72  01  -0.X47962  Cl  -0.926632  00  -0.710412  00  -0.391412-00  -0.231612-00  -0.132022-00  -0.772082-01  -0.t6ll6£-0l 
0.308672-01  0.733992-01  0.111972-00  0.146712-00  0.179192-00  0.210442-00  0.241102-00  0.27IS62-00  0.302002-00 

0.3'X2402-00  0.362332-00  0.392022-00  0.420132-00  0.446032-00  0.468612-00  0.466312-00  0.496362-00  0.302692  00 

•.498972-00  0.486032-00  0.469422-00 

-0*374332  04  -0.337012  04  -0.173922  04  -0.346732  03  -0.160362  03  -0.342732  02  -0.217092  02  -0*101092  02  -0.336082  01 

-0.321902  01  -0« 212942  01  -0.133222  Ot  -0.129962  01 

-0.374332  04  -0.337022  04  -0.173942  04  -0.346832  03  -0.160412  03  -0.342982  02  -0.217232  02  -0.I0SI62  02  -0.336692  Ot 

-0.322362  vt  -0.213302  01  -0.133SX  01  -0.130262  01  -0.962232  00  -0.813962  00  -0.714332  00  -0.641382  00  -0.386872  00 

-0.343142  00  -0.312472  00  -0.486412-QO  -0.466172-00  -0.430022-00  -0.436602-00  -0.423692-00  -0.416082-00  -0.40749C-00 
-0.399302-00  -0.391762-00  -0.3al94£-Q0  -0.375742-00  -0.366912-00  -0.337222-00  -0.346332-00  -0.334782-00  -0.322032-00 
-0.308322-00  -0.294322-00  -0.263392-00 

0.17X992  04  0.134842  04  0.808X36  03  0.231012  03  0.734372  02  0.247072  02  0.974102  01  0.440922  01  0.223302  01 

•.12403C  01  0.738196  00  0.462222-00  0.334322-00 

0.171^92  04  0.134842  04  0.808222  03  0.2'JI062  03  0.734802  02  0.247192  02  0.974732  01  0.441302  01  0.223732  Ot 

0.t2420C  01  0.739462  00  0.463I62-00  0.333132-00  0.193722-00  0*125842-00  0.760482-01  0.386342-01  0.906832-02 

-8.tS434C-0l  -0.367S9C-0I  -0.360932-01  -0.733182-01  -0.698242-01  -0.103322-00  -0.120932-00  -0.136262-00  -0.131382-00 
-••166892-00  -0.162082-00  -0.196942-00  -0.211142-00  -0.224232-00  -0.233642-00  -0.244722-00  -0.230762-00  -0.233112-00 
-0.231222-00  -0.244792-00  -0.236472-00 

0.374312-03  -0.948732  03  -0.148312  04  -0.916372  03  -0.441452  03  -0.217712  03  -0.II392C  03  -0.664982  02  -0.406202  02 
-0.26X222  92  -0.173192  02  -0.121632  02  -0.989892  Ot 

-0.374912-09  0.948732  03  0.148312  64  0.916372  03  0.441452  03  0.217712  03  0.113922  03  0.664982  02  0.406202  02 

0.261222  02  0.173192  02  0.171632  02  0.989692  01  0.633392  01  0.474042  01  0.339132  Ot  0.273412  01  0.213042  01 

••169602  01  0.I28702  Ot  0.996292  00  0.770872  00  0.390932  00  0.443312-00  0.319212-00  0.212702-00  0.119732-00 

••37370C-FX  -0.333332-01  -0.100622-00  -0.136072-00  -0.207792-00  -0*249262-00  -0. 281712-00  -0.304212-00  -0.313922-00 
-••3IO31‘-‘00  -0.30S39C-Q0  -0.289062-00 

••  0.440142  03  0.680032  03  0.423132  03  0.204772  03  0.100972  03  0.337482  02  0.308232  02  O.I6S20S  02 

0.120972  02  0.210992  01  0.362932  Cl  0.438172  01 

-••  -0.440142  03  -0.688032  03  -0.423102  03  -0.204772  03  -0.100972  03  -0*337462  02  -0*308232  02  -0.188202  02 

-0.120972  02  -0.010992  0|  -0.962932  01  -0.438172  01  -0.294342  01  -0.220212  01  -0.167612  01  -0.129372  Ot  -0.101372  Ot 

-•••06372  00  -0.64C222  00  -0.327962  00  -0.440302-00  -0.376002-00  -0.329262-00  -0.296362-00  -0.274922-00  -0.263392-00 

-••260032-00  -0.266692-00  -0.260472-30  -0.302002-00  -0.330072-00  -0.3663^2-00  -0.406032-00  -0.434092-00  -0.302762  00 
•••9SI032  00  -0.397902  00  -0.630762  00 
••.7030 

0.29X732  02  0.343602  01  0.128902  01  0.618922  00  0.293662-00  0.964912-01  -0.315432-01  -0.1033S2-00 

-••134662  02  -0.153362  01  -0.943132  00  -0.226872-00  -0.640292-01  -0.109122-01  0.293672-01  0.306682-01 

THCTASVtO  *  3 

lMrCA«l20IAr2  iMFOAMATtON  OF  VHL.VHT.get.UST. VBL.VOt  AT  P  2 

0.  0.349072-02  0.872682-02  0.137082-01  0.244352-01  0.34907C-01  0.471232-01  0.610872-01  0.767962-01 

••942492-01  0.II349C-00  0.121442-00  0.197082-00  0.181322-00  0.207702-00  0.235622-00  0.265292-00  0.296712-00 

0.329872-00  0.364702-00  0.399692-00  0.434592-00  0.469502-00  0.5044X2  00  0.539312  00  0.574222  00  0.609I3C  00 

0. 644042  00  0.670942  00  0.713652  00  0.740762  00  0.783662  00  0.016572  00  0.653482  00  0.888382  00  0.923292  00 

0.958202  00  0.98147C  00 

0.191532  04  0.172312  04  0.109092  04  0.493772  03  0.199242  03  0.835942  02  0  384752  02  0.196112  02  0.110282  02 

0.679162  01  0.454052  01  0.395762  01  0.249342  01  0.200922  01  0.168702  01  0.146402  01  O.I306IC  01  0.118932  01 

0.1I0II2  01  0.103302  01  0.981662  00  0.941772  00  0.909932  00  0.883642  00  0.061872  00  0.842732  00  0.825482  00 

0.809222  00  0.793212  00  0.776792  00  0.759362  00  0.740422  00  0.719592  90  0.696692  00  0.671742  00  0.64307E  00 

0.617252  00  0.398462  00 

-0.071322  03  -0.703832  93  -0.4961*02  03  -0.224292  03  -0.902312  02  -0.375882  02  -0.170462  02  -0.e4533e  01  -0.454A76  01 

-0.261252  Ot  -0.158392  01  -0.131702  01  -0.643962  00  **0.417072-00  -0.263662-00  -0.154972-00  -0.735902-01  -0.966062-02 

0.431152-01  0.888492-01  0.128472-00  0.164412-00  0.198162-00  0.230696-00  0.262562-00  0.294132-00  0.323442-00 

0.356342-00  0.366402-00  0.415232-00  0.441742-00  0.464892-00  0.483412-00  0.495692-00  0.500992  00  0.497602-00 

0.485072-00  0.471622-00 

-0.952832  03  -0.857156  03  -0.542632  03  -0.245572  C3  -0.990732  02  -0.415602  02  -0.191232  02  -0.974572  01  -0.347692  01 

-0.337312  01  -0.225432  01  -0.196462  01 

-0.952832  03  -0.857252  03  -0.542792  03  -0.245/06  03  -0.991582  02  -0.416112  02  -0.191562  0?  -0.976652  Ot  -0.549J62  01 

-0.338422  01  -C.-F26322  01  -0.197302  01  -0.124462  01  -0.100272  01  -0.642312  00  -0.731726  03  -0.652/72  00  -0.59472E  00 

-0.550922  00  -0.517172  00  -0.491762-00  -0.472072-00  -0.456^62-00  -0.443572-00  -0.432802-00  -0.423*52-00  -0*416032-00 

-0.407102-00  -0.399292-00  -0.391262-00  -0-382712-00  -0*373382-00  -0*303102-00  -0*351752-00  -0.339362-00  -0*326002-00 

-0.3I22IE-00  -0.302826-00 

0.433462  03  0.389912  03  0*246762  03  0*111552  03  0.448682  02  0.166062  02  0.647302  01  0.420192  01  0.225792  01 

0.129752  01  0.766412  00  0.653792  00 

0.433462  03  0.369962  03  0.246832  03  0.111612  03  0.449062  02  0.187102  02  0.848702  01  0.421092  01  0.226392  01 

0.130162  01  0*769532  00  0*656572  00  0*321242-00  0.208142-00  0.131742-00  0.774092-01  0.367792-CI  0.463062-02 

-0.215722-01  -0.444822-01  -0.643592-01  -0.024112-01  -0.993902-01  -0.115772-00  -0.13166F-00  -0.147796-00  -0.163622-00 

-0*179272-00  -0.194552-00  -0*209142-00  -0*222632-00  -0.234442-00  -0*243922-00  -0.250372-00  -0.253102-00  -0.251542-00 

-0.245352-00  -0.238642-00 

0.370222-05  -0.236962  03  -0*374142  03  -0*302612  03  -0*187472  03  -0.109592  03  -0.653562  02  -0.405792  02  -0.262692  02 

-0.176602  02  -0.122656  02  -0.107172  02 

-0.370222-05  0.23696C  03  0.374142  03  0.302812  03  0.187472  03  0.1099^  03  0.653562  02  0.405792  02  0.262692  02 

0.176602  02  0.122652  02  0.107172  02  0.640132  01  0.4V700E  VI  0.361092  01  0.276702  01  0.213922  01  0.166222  01 

0.129252  01  O.lOOOti  01  0.774292  00  0.594222  00  0*446562-00  0*322552-00  0.216162-00  0.I233DE-P0  0.412112-01 

-0.319022-01  -0.970812-01  -0.154752-00  -0.204842-00  -0*246842-00  -0.279962-00  -0.303306-00  -0.315946-C'  -0 .31 728JP-00 

-0.307242-00  -0*294512-00 

0.  0*109932  03  0*173572  03  0*140472  03  0*869492  02  0.506172  02  0.302962  02  0.168032  02  0.121662  02 

0.817562  01  0.567702  01  0*496042  01 

-0.  -0.109932  03  -0.173372  03  -0*140472  03  -0*669492  02  -0.500172  0?  -0.302962  02  -0.188032  02  -0.121662  02 

-0.817562  01  -0.567702  01  -0.496042  01  -0*296632  01  -0*221552  01  -0.168462  Ot  -0.130102  01  -0.101692  01  -0.606232  00 

-0.049162  00  -0*528262  00  -0*440142-00  -0*375462-00  -0*326422-00  -0.295222-00  -0.273472-00  -0.261622-00  -0.258712-00 

-0*264172-00  -0.2776^-00  -0.298816-00  -0*327382-00  -0.362832-00  -ft. 404302-00  -0.450482-00  -0.499522-00  -0.549042  00 

-0.896242  00  -0*624932  00 
••*7100 

0.818692  02  0.410842  01  0*144772  01  0*716922  00  0.375262-00  0.169122-00  0*336542-01  -0.42053C-0I 

-0.100272  02  -0.16S23E  01  -0.609186  00  -0.265172-00  -0.II6I26-00  -0.410252-01  0.36336E-0J  0.221822-01 

(similar  print  out  for  R  *  0.715,  0.721 .  0.9999) 
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HUS-XCmC 

rhcrA$u0O  3 

iNTCMMeOIAre  IK^QPMATION  of  PM|«U»'H«VHH»UHL«UHr.VHU.VHl.U01.«UOT.VHL«Var  AT  ois 


e*  o«i745ac*04p  o«S23eie>o2  o.iOAr^cf.oi 

0»7eSAte-0t  0.9S994£~01  0*I1S19C'>00  0.tJM4e-0d 

0*29S4ee*00  0«33l62e-00  0.366S2C-00  0.40143F-00 

0«610S7E  00  0«64S7eC  00  0«6a064e  00  0.715:»9£  00 

C*92904e  00  0*9S994E  00  0.994e5C  00  0.99973C  30 

-0«7t780E  04  ''0.64636I:  04  «0.3370eC  04  -O.I0424C  04 
-0«l7798e  01  «0«27934e>00  0«4|430C>00  0*7J744k'  00 

0«90ia6C  00  0«87779e  00  0.8S430C  00  0«83312C  00 

0*73e27e  00  0«7237e£  00  0.70843C  CO  0..69180e  00 

0«9S282e  00  0«S24ee6  00  0.49737E-00  0.493ft2e-00 

e«33626C  04  0.302eoe  04  0«1S83S£  04  0.4a87SC  03 

0«l3l99e  01  0.61706e  00  0.29362E*00  0.t4470£-00 

0«tOS83e-00  0«t2a92C-00  0.tS383C*-00  0.l79l6t:*00 

0.342Ste-00  0,371ll£-00  0.39926C-00  0.42623e-00 

O.S0I07E  00  o.4ee26£-oo  o.4a6a9C-oo  o.4a32ie‘-oo 

0*361 1  IE  04  C*329ia£  04  0*l7004e  04  0.S2446E  03 

0*8999SE  00  O.I4064E-CO  -0*20a67F*00  <-0*37|42£-^0 

0«36llte  04  0*32S17£  04  0*t7002E  04  0.92437E  03 

0«a9477E  00  0*14042E*00  -0«20827£-00  *0* J70t>7&-00 

>>0*492S7£>00  -0.44036E-00  >0*42a46E-00  ^0*4l77t£-00 
-0*369SOe-00  -0*362t4E-00  *0*3E436E'*00  *0*34594£*00 
•0*27999E>00  -0*26t92£*00  -0*24fll2£«00  -0*2a624E>*00 
•>0*t69t6E  04  -0*I9233E  04  -0*7966eE  03  -0.24S90E  03 
•0*66492e  00  -0*31009C-00  -0*t47e6E«00  •0«7267e£-0t 
>'0*I69I6E  04  -0*13233£  04  -0.79661E  03  •0*249tf6E  03 
-0*663I4E  00  -0«3t019E-00  -0«I4797E-00  -0*72712E-0t 
-0*93tlOE-Ot  ->0.6467ae-0l  -0*77149E-01  -0.a9827£-0t 
<-0*17143E-00  -0*18969e-00  -0*l997lE>ao  -0*2t314e*00 
*0*290126-00  *0*24365£*00  *0*23290E*00  *0«23107C*00 
0.38207E-OS  *0*99016E  03  *0*1493IE  04  -0*926676  03 

-0.26926E  02  *0*I7791£  02  -0*t2392e  02  -0*ae2J9E  01 

-0*382876-’0S  0*99016E  03  0*|4031E  04  0*926676  03 

0.26S26E  02  0*177916  02  0.:2392£  02  0*a0239C  01 

0*160086  01  0*13068e  01  0*10I07C  01  0*701/76  00 

0.379S2E-01  *0«398466*01  -0*101366-00  -0*198996-00 
-0»31406E-C0  -0*3024^-00  -0«20Q37£-00  -0*276946-00 
0*  0*440806  03  0«69268£  03  0*429886  03 

0.12283E  02  0*e239aC  01  0*S7168£  Ot  0*400906  01 

-0*  -0*440806  03  -0*692686  03  -0.4290AC  03 

-0*122a3E  02  -0*623986  01  -0*571686  01  -0*408906  01 

-0*8I09IE  00  -0*698306  00  -0*936196  CO  -0*447186-00 
-0*264796-00  -0*270476-00  -0«2a421£-00  -0*30562£-0d 
-0*9StOOC  00  -0*99970£  00  -0*63470£  00  -0«639bSC  00 
8**6930 

0*7S944E  01  0*27396£  01  0*132246  01  0*77237£  00 

-0«>3?.$7ie  01  -0*|2S13F  01  -0*978606  00  -0*312296-00 
THETASuaO  3 

INTEPFEOIATE  thFORMATION  OF  PH  f  *Ut-H*  VHi«*  UHt.  •WHT  •  VHL 
0*  0*349076-02  0*872686-02  0«19708C-0I 

0*942496-01  0*113456-00  0*134396-00  0*145796-00 

0*329876-00  0*3647A6-00  0*39969e-00  0*434996-00 

0*644046  00  0*67894C  00  0*713896  00  0*748/66  00 

0*990206  00  0*993116  00  0*100976  01 

-0*179386  04  -0*15796£  04  -0*100266  04  -0*4SI69£  03 

-0*260826  01  -0*936656  00  -0*844736-01  0«|7869£-00 

0*807036  00  0«79766C  00  0*789996  00  0*773*>9E  00 

0*696646  00  0*681626  00  0*669096  00  0«64673E  00 

0*497416-00  0*470226-00  0*460776-00 

0*829886  03  0*747496  03  0*474626  03  0«2l4ll6  03 

0*172136  01  0*931496  00  0*53022C  OC  0*407016-00 

0*195976-00  0«1739?£-00  0«19409€-00  0*216396-00 

0*373966-00  0*401036-00  0*427466-00  0*451906-00 

0*489456-00  0*468666-00  0,459256-00 

0*885666  03  0*797196  03  0*906386  03  0*228176  03 

0*131936  01  0*473956-00  0,427596-01  -0*904706-01 

0*885686  03  0*797656  03  0*906236  03  0*228096  03 

0*131906  01  0*472086-00  0*429606-01  -0«900I2C-0I 

-0*409226-00  -0*4002n:-00  -0*394176-00  -0*387696-00 
-0*347866-00  -0*340166-00  -0«33170£-C0  -0*3?234C-00 
-0*247016-00  -0*233366-00  -0*228626-00 
-0*419106  03  -0*377916  03  -0*239736  03  -0*108166  03 
-0*870686  00  -0*471346-00  -0*268396-00  -0*206066-00 
-0*419106  03  -0*17'’47E  03  -0*239696  03  -0*10ai0£  03 
-0*867836  00  -0*469466-00  -0*267146-00  -0«?09036-00 
-0*781156-01  -0.870706-01  -0*973136-01  -C  #  I08456-.00 
-0*186546-00  -0*200136-00  -0*213196-00  -0*225246-00 
-0*243056-00  -0*232596-00  -0*227e<6-00 
0*387256-09  -0.237646  03  -0.377996  03  -0*108956  03 
-0*182046  02  -0*126496  02  -0*903016  01  -0*767146  01 
-0*307256-05  0.23764F  03  0*377996  03  0*308556  03 

0*182046  02  0*126496  02  0*903016  01  0«767|46  01 

0*133106  01  0*102936  01  0*796366  00  0*610636  00 

-0*324676-01  -0*985366-01  -0*196606-00  -0*206566-00 
-0*301466-00  -0*279686-00  -0*269976-00 
0*  0*110256  03  0*175396  03  0*143136  03 

0*842796  01  0*989326  01  0*418066  01  0*395316  01 

-0*  -0*119256  03  -0*175356  03  -0,143l3€  03 

-0*842796  01  -0*9bS32£  01  -0*418066  01  -0*355316  01 
-0*669506  00  -0*544846  00  -0*493996-00  -0*387336-00 
-0*271816-00  -0*289106-00  -0*306066-00  -0*334256-00 
-0*992026  00  -0«C304S6  00  -0*642466  00 
Mt^6900 

04903376  01  0*234736  01  0*118406  01  0*600566  00 

-0*232396  01  -0* 107806  01  -0*523086  00  -0*276796-00 

(Similar  print  out  for  R»  0.685,  0.680,  0.2001) 


0*174546-01  0*26ia6£-0l  0*36652C-01  0*468706-01 

0*139796-00  0*183266-00  0*209446-00  0*237376-00 

0*436346-00  0*471256-00  0.506I5C  00  0*541066  00 

0*750506  00  0*785416  00  0*820326  00  0*855226  00 

-0*297566  03  -0.949I9C  02  -0.340526  02  -0*131716  02 
0*770666  00  0*939526  00  0*952566  00  0*943856  00 

0*814136  00  0*796986  00  0*761276  00  0*766566  00 

0*673516  00  0*653276  00  0*630926  00  0*606506  00 

0*139876  03  0*449516  01  0*164396  02  0*665576  01 

0*179626-00  0*577046-01  0*570886-01  C*679pl6-01 

0*204946-00  0*231276-00  0*258216-00  0.28580F-00 

0*451106-00  0*472756-00  0*489876-00  0*501106  00 

0*149726  03  0*477626  02  0*171366  02  0*662886  01 

-o*3eai6£-oo 

0.14967C  03  0*477416  02  0*l7126fc  02  0*662336  01 

-0*387256-00  -0*471936-00  -0*478376-00  -0*473896-00 
-0*408076-00  -0*399366-00  -0*391376-00  -0*383896-00 
-0*336696-00  -0*326466-00  -0*315216-00  -0*302926-00 

-0*703766  02  -0*228196  02  -0*827246  01  -0*334976  01 
-0*652056-01 

-0*703556  02  -0*22609h  02  -0*826746  01  -0*33«706  01 
-0*651326-01  -0*209856-01  -0*286696-01  -0*341216-01 
-0.I0272F-00  -0*tl5A9F-00  -0*12935^-00  -0*143136-00 
-0*225516-00  -0*23626£-00  -0*2«474c-00  -0*250266-00 

-0*447536  03  -0*220916  03  -0*117676  03  -0*675196  02 
-0*836386  01 

0*447536  03  0*220916  03  0*1176/6  03  0*675196  02 

0*836386  01  0*48l38b  Ot  0*364666  01  0*279836  01 

0*599056  00  0*449l6L-00  0*J2322b-00  0*215216-00 

-0*208646-00  -0*249616-00  -0*201736-00  -0*303546-00 

0*207596  03  0*102466  03  0*545626  02  0*312966  02 

0*387256  01 

-0*207596  03  -0«10248F  03  -0»54S62£  02  -0*312966  02 

*-0*J8725F  01  -0*22364t  01  -0*170226  01  -0*131596  Ot 

-0*381896-00  -0*334436-00  -0*301026-00  -8*279226-00 
-0*334306-00  -0*369616-00  -0*410816-00  -0*455676-00 


0*480596-00  0*313006-00  0*19714fc-00  0*131136-00 

-0*183016-00  -0*115816-00  -0*787706-01  -0*593636-01 


VHT*Ultt.*U8T«  vei.«yor  at  P16 


0*244356-01 

0*349076-01 

0*471256-01 

0*610876-01 

0*181526-00 

0*207706-00 

0.235626*00 

0*265296-00 

0*469506-00 

0*504416  00 

0*539316  00 

0*574226  00 

0*783666  00 

0*818576  00 

0.85346F  00 

0*888366  GO 

-0*178256  03 

-0*712526  02 

-0*300966  02 

-0*133376  02 

0*595516  00 

0*717526  00 

0*777826  00 

0*603036  OQ 

0*761016  00 

0*748566  00 

0./36146  00 

0*723556  00 

0*626306  00 

0*603736  00 

0*57Vt0F  00 

0*552746  00 

0*647996  02 

0*341926  02 

0*147256  02 

0*679696  01 

0*215266-00 

0*  162626-00 

0*140506-00 

0*  136096-00 

0*240106-00 

0*265196-00 

0*291276-00 

0*318226-00 

0*473236-00 

0*490166-00 

0*501366  00 

0*505576  00 

0*900576  02 

0.36005E  02 

0*152116  02 

0*674266  Cl 

0*899816  02 

0*359616  02 

0*151866  02 

0*672936  01 

-0*299786-00 

-0*361046-00 

-0*391 196-00 

-0*404066-00 

-0*381106-00 

-0.37471E-00 

-0*368276-00 

-0*361 746-00 

-0*311976-00 

-0*300546-00 

-0*288106-00 

-0*274«.>e-00 

-0*428426  02 

-0.172786  02 

-0*744236  01 

-0*343616  01 

-0*428066  02 

-0*l7257t  02 

-0*743016  OS 

-0*342686  01 

-0*108366-00 

-0.618296-01 

-0.706646-01 

-0*6841 lF-01 

-0*120306-00 

-0*132756-00 

-0*l457lF-00 

-0*159096-00 

-0*235726-00 

-0*244016-00 

-0.249436-00 

-0.251376-00 

-0*192096  03 

-0*112636  03 

-0*672696  02 

-0*418006  02 

0*192096  03 

0*  1 12636  03 

0*672696  02 

0*418006  02 

0*491796  01 

0*372246  01 

0*2852lk  01 

0*220456  01 

0*45045e.00 

0*330726-00 

0*221316-00 

0*126026-00 

-0*247976-00 

-0*280076-00 

-0*302036-00 

-0*3131 16-00 

0*890936  02 

0*522256  02 

0*311836  02 

0*193696  02 

-0.890936  02 

-0*522256  02 

-0.3ll83>i  02 

-0*I9369C  02 

-0*226476  01 

-0*173726  01 

-0*134176  01 

-0*105006  01 

-0*338016-00 

-0*304576-00 

-0*202096-00 

-0*269766-00 

-0*369056-00 

-0.409486-00 

-0.454122-00 

-0*S0107F  00 

0*411966-00 

0*243966-00 

0*132776-00 

0*693736-01 

-0.152316-00 

-0*865126-01 

-0.500426-01 

-0*309036-01 

O.82033C-O1 
0*267046-00 
0.57S976  00 
0*890136  00 

-8.514496  01 
0*924946  00 
0*752386  00 
0*580296  00 

0*289536  01 
0.050906-01 
0*313966-00 
0. 5051 76  00 

0*258976  01 

0*250696  01 
-0*484276-00 
-0-3/6886-00 
-0.289756-00 

-0*145746  01 

-9*145586  Ot 
-0.4271 lE-Ot 
-0.l5/l8r-08 
-0.252246-00 

-0*4124*6  02 

0*412476  02 
0.21&276  01 
0.I2I05E-0U 
-0*31«<6C-00 

0*1911 16  02 

-0*191116  02 
-0*103156  01 
-0,287466-00 
-0*503476  00 


0*767966-01 
0.2787IC-00 
0*609136  00 
0*923291;  00 

-0*60167e  01 
0*810716  00 
0*710496  00 
0.525256  00 

0*333356  01 
0*142406-00 
0*345766-00 
0*501766  00 

0*304266  01 

0*303446  01 
-0*407306-00 
-0*355006-00 
-0.260996-00 

-0*168576  01 

-0*168126  01 
-0*715426-01 
-0.172776-00 
-0*249336-00 

-0*270726  02 

U.270726  02 
0*171246  01 
0*420386-01 
-0*312886-00 

0*125396  02 

-0*125396  02 
-0.833536  00 
-0*266556-00 
-0*547956  00 
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TASUBO  3 

CRXCOtATC  INFORNATICN  OF  PH I • UKH • VHh*UHL •UHT » VHL • 
't745AC'02  O.S2361C--02  O.10472E-0I  0.174546-01 

95094E-01  0«ll519e>00  0.I3374£*00  0.t58e3£-00 

'33I62C-0C  0.36652E-00  0.40t43E-00  0.43634k-00 

I64578E  00  0.0e069£  00  0.71SS9E  00  0.7S0S0E  00 

I9S994E  00  0.99370E  00 

'S6907E  00  0.20416E-00  0.1I276E-00  0.738S9E-01 

14S13E-01  0.12I06E-01  0.10435E-0t  0.67917£«02 

1422146-02  0.382I4E-02  0.34910E-02  0.32137E-02 

'21810E-02  0«20710£-02  0.197I6E-02  0.18619E-02 

I14812E-02  0*  14  326£-'02 

i6eS03E-02  •0.3l023£-0t  -0.380486-01  -0.3057SE-01 
1669756-02  -0.55895E-02  -0.48172E-02  -0.40560E-02 
1192186-02  -0.173266-02  -0.157566-02  -0.144326-02 
940916-03  -0.885936-03  -0.63095S-03  -0.790276-03 
579936-03  -0.553366-03 

285216-00  -0.102326-00  -0.56514E-01  -0.370176-01 
727346-02  -0.606726-02  -0.522966-02 
1265226-00  -0.10232E-00  -0.565156-01  -0.370186-01 
727396-02  -0.606786-02  -0.523016-02  -0.440666-02 
211606-02  -0.191556-02  -0.174996-02  -0.161096-02 
I1093A6-02  -0.103826-02  -0.988516-03  -0.94346E-33 
742626-03  -0.718276-03 

343336-02  0.155486-01  0.190696-01  0.153246-01 

335666-02  0.280136-02  0.241426-02 

343336-02  0.155486-01  0.190706-01  0.153246-01 

335696-02  0.280156-02  0.241456-02  0.203296-02 

963326-03  0.868476-03  0.789796-03  0.723416-03 

471696-03  0.444136-03  0.419086-03  0.396166-03 

290766-03  0.277446-01 

140346  03  -0.376006  02  -0.127456  02  -0.505626  01 
177806-00  -0. 123576-00  -0.916886-01 
140346  03  0.376606  02  0.127456  02  0.505626  01 

177866-00  0.123576-00  0*916886-01  0*649996-01 

148386-01  0.121286-01  0.100936-01  0*852646-02 

383406-02  0.343996-02  0.310206-02  0*281006-02 

167266-02  0. 1 55406-02 
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390406-03  -0.334796-03 
6CPAT60  VA»VT  IN  Ml C-ZONE 

107966-01  0.230586-02  0.949466-03  0.491276-03 

496946-02  -0.103636-02  -0.412246-03  -0.206086-03 
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Figure  17  (Continued) 


53 


APPENDIX  D 

FORTRAN  LISTING  OF  COMPUTER  PROGRAM 

The  FORTRAN  listing  of  the  computer  program  is  included  in  the  follow¬ 
ing  pages^.  The  computer  also  uses  the  subroutine  MATINV  and  the  function 
SIMPUN. 
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